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ABSTRACT 

We present SED modeling of a sample of 34 K-selected galaxies at z ^ 2.3. These galaxies have 
NIR spectroscopy that samples the rest-frame Balmer/4000A break as well as deep photometry in 
thirteen broadband filters. New to our analysis is IRAC data that extend the SEDs into the rest-frame 
NIR. Comparing parameters determined from SED fits with and without the IRAC data we find that 
the IRAC photometry significantly improves the confidence intervals of r, Ay, Mgtar, and SFR for 
individual galaxies, but does not systematically alter the mean parameters of the sample. We use 
the IRAC data to assess how well current stellar population synthesis codes describe the rest-frame 
NIR SEDs of young galaxies where discrepancies between treatments of the TP-AGB phase of stellar 
evolution are most pronounced. The models of Bruzual & Chariot (2003), Maraston (2005), and 
Chariot & Bruzual (2008) all successfully reproduce the SEDs of our galaxies with < 5% differences 
in the quality of fit; however, the best-fit masses from each code differ systematically by as much 
as a factor of 1.5, and other parameters vary more, up to factors of 2-3. A comparison of best-fit 
stellar population parameters from different SPS codes, dust laws, and metallicities shows that the 
choice of SPS code is the largest systematic uncertainty in most parameters, and that systematic 
uncertainties are typically larger than the formal random uncertainties. The SED fitting confirms 
our previous result that galaxies with strongly suppressed star formation account for '-^ 50% of the 
K-bright population at z ~ 2.3; however, the uncertainty in this fraction is large due to systematic 
differences in the SSFRs derived from the three SPS models. 

Subject headings: infrared: galaxies — galaxies: fundamental parameters — galaxies: evolution — 
galaxies: stellar content — galaxies: high-redshift 
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1. INTRODUCTION 

The observed flux from a galaxy as a function of 
wavelength or frequency, its spectral energy distribu- 
tion (SED), represents the integrated light of its stel- 
lar populations and is a valuable tool for determining 
properties such as its age, and current mass contained 
in stars. Interpreting observations of galaxy SEDs re- 
quires creating synthetic spectra using our best under- 
standing of stellar spectra, stellar evolution, and dust 
absorption. This process, known as stellar population 
synthesis (SPS) was first implemented by Tinsley (1967, 
1972) and has subsequently been refined into increas- 
ingly sophisticated models (e.g., Bruzual 1983; Renzini 
& Buzzoni 1986; Bruzual & Chariot 1993; Worthey 1994; 
Fioc & Rocca-Volmerange 1997; Leitherer et al. 1999; 
Bruzual & Chariot 2003; Maraston 2005; Bruzual 2007; 
Conroy et al. 2008). 

With the advent of deep multiwavelength extragalactic 
surveys and the refinement of photometric redshift tech- 
niques, large, statistically complete samples of galaxies 
with broadband SEDs are now available for studying the 
evolution of galaxy properties in a systematic way up to 
z ~ 2-3 (e.g., Forster Schreiber et al. 2004; Labbe et al. 
2005; van Dokkum et al. 2006; Papovich et al. 2006; 

^ Department of Astronomy, Yale University, New Haven, CT, 
06520-8101; adam.muzzin@yale.edu 

^ Hubble Fellow, Carnegie Observatories, 813 Santa Barbara 
Street, Pasadena, CA, 91101 

^ Department of Astropliysical Sciences, Princeton University, 
Princeton, NJ, 08544 

Leiden Observatory, Leiden University, PO Box 9513, 2300 RA 
Leiden, Netherlands 



Fontana et al. 2006; Daddi et al. 2007; Wuyts et al. 
2007; Perez- Gonzalez et al. 2008; Drory & Alvarez 2008; 
Marchesini et al. 2008, and many others). The demand 
for SPS codes that provide accurate synthetic SEDs for 
a given parameter set is now greater than ever as they 
are the tool for extracting astrophysical information from 
these data. 

Using the current generation of SPS models these stud- 
ies have shown that a large fraction of the stellar mass 
(Mstar) in massive galaxies is assembled by z ~ 2 (e.g., 
Fontana et al. 2006; Perez- Gonzalez et al. 2008; March- 
esini et al. 2008), and that a significant fraction of these 
galaxies are already passively evolving (e.g., Labbe et al. 
2005; Papovich et al. 2006; Kriek et al. 2006). Unfortu- 
nately, the current SPS models do have significant chal- 
lenges to overcome (see e.g., Conroy et al. 2008), making 
it unclear how robust some of these results are (see e.g., 
the discussion in Marchesini et al. 2008). Furthermore, 
it is well known that stellar population parameters de- 
rived from fitting broadband SEDs suffer from system- 
atic uncertainties caused by the fact that metallicity, the 
initial mass function (IMF), and dust extinction law can- 
not be constrained by broadband SEDs alone. Instead, 
these parameters must be assumed a priori when creat- 
ing the set of synthetic SEDs used in the fitting. More 
recently it has become clear that the treatment of SPS 
itself (e.g., the choice of isochrones, stellar libraries, and 
integration method) can result in synthetic spectra that 
look strikingly different for a similar set of input parame- 
ters (e.g., Maraston 2005, Maraston et al. 2006; Bruzual 
2007; Conroy et al. 2008) and that the choice of SPS 
method is now an additional systematic error in the in- 
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terpretation of galaxy SEDs. 

One of the major challenges for SPS remains to be a 
complete treatment of the thermally-pulsating asymp- 
totic giant branch (TP-AGB) stars. This population 
of stars are notoriously challenging to model because 
they have complex physics such as thermal pulses, evolv- 
ing chemical compositions, rapid mass-loss rates, and 
self-obscuring dust shells (see e.g., Marigo ct al. 2008, 
and references therein). Furthermore, their lifetimes are 
short which means they are rarely found in clusters and 
therefore the majority of empirical spectra come from 
field stars of unknown metallicity (e.g., Langon & Wood 
2000, Langon & Mouchine 2002). Nevertheless, they can- 
not be ignored in SPS because they are bright in the rest- 
frame NIR and can dominate the total NIR luminosity 
of young stellar populations with ages between 0.2 to 2.0 
Gyr (e.g., Maraston 2005, Bruzual 2007). The treatment 
of TP-AGB stars in SPS modeling is now a key issue in 
the study of stellar masses, because the rest-frame NIR 
also traces a galaxy's old stellar population and therefore 
represents its integrated star formation history (SFH) . If 
galaxies undergo multiple bursts of star formation (SF) 
at widely spaced epochs, the rest-frame NIR becomes the 
critical wavelength range for determining stellar masses. 
For young galaxies at high rcdshift, SPS codes with dif- 
ferent treatments of TP-AGB evolution can produce stel- 
lar masses and ages that systematically differ by roughly 
a factor of 2 for identical SEDs (e.g., Maraston et al. 
2006; Kannappan & Gawiser 2007; Wuyts et al. 2007; 
Bruzual 2007). Indeed, based on discrepancies between 
kinematically-derived, and photometrically-derived stel- 
lar masses of galaxies a.t z ^ 1, van der Wei et al. (2006) 
have even argued that the stellar masses of young, early- 
type galaxies are better determined by SED fits that ex- 
clude the rest-frame NIR. 

Recent papers by Maraston et al. (2006) and Bruzual 

(2007) have suggested that Spitzer observations could 
be important for resolving some of the differences be- 
tween SPS codes, particularly the TP-AGB treatment. 
IRAC observes at 3.6/im, 4.5/im, 5.8/im and 8.0/im, ef- 
fectively the rest- frame NIR of galaxies at z ~ 2. At this 
redshift the universe is young, ~ 3 Gyr, and therefore 
the majority of observed galaxies contain the young stel- 
lar populations that should be dominated by TP-AGB 
stars; however, they have not yet acquired an underlying 
old stellar population that must be accounted for when 
studying TP-AGB effects in the rest-frame NIR in low 
redshift galaxies (e.g., Riffel et al. 2008). Using a sam- 
ple of 7 galaxies at lA < z < 2.7 from the GOODS 
survey, Maraston et al. (2006) have argued that their 
models produce better fits to young galaxies than the 
Bruzual & Chariot (2003, hereafter BC03) models and 
this is caused by their empirically calibrated fuel con- 
sumption treatment of the TP-AGB phase. Using the 
next generation of the BC03 code that includes improved 
evolutionary tracks for TP-AGB stars from Marigo et al. 

(2008) , Bruzual (2007) have fit the same 7 galaxies and 
showed that they now obtain parameters closer to those 
in Maraston et al. (2006). 

In Kriek et al. (2006; 2007; 2008a, hereafter K08) we 
presented SED modeling of a sample of 36 K-selected 
galaxies at z ~ 2.3 that was based on extensive NIR 
spectroscopy covering the JHK bands obtained from 
GNIRS on Gemini, as well as broadband photometry. 



These data constrain the galaxy's SEDs from the rest- 
frame UV to optical and K08 compared fits with and 
without the NIR spectroscopy to determine how pa- 
rameters determined from fitting broadband data alone 
compare with those determined with the higher reso- 
lution spectroscopy and spectroscopic redshifts. Re- 
cently, we have obtained IRAC observations of this sam- 
ple and this allows us to extend the SED modeling 
into the rest-frame NIR. With deep photometry in 13 
broadband fihers (UBVRIz' JHK+IRAC) as well as high- 
resolution NIR spectroscopy that covers the rest-frame 
Balmer/4000A break, this sample of galaxies are likely to 
have the best-constrained SEDs of young massive galax- 
ies at z ^ 2.3 mitil more powerful optical and mid- 
infrared spectroscopic capabilities become available from 
ground-based 20-30m telescopes and JWST, respectively. 

In this paper we explore how well the properties of 
these galaxies can be derived from their SEDs, with a 
focus on the rest-frame NIR. Specifically, we concentrate 
on four issues, 1) Improvement in Constraints and Sys- 
tematic Errors: Does IRAC data improve the constraints 
on stellar population parameters for z ~ 2.3 galaxies, 
or does it increase the uncertainties because the mod- 
els fail to correctly reproduce the rest-frame NIR SED? 
Are there systematic errors in the stellar population pa- 
rameters determined without the rest-frame NIR data? 
2) Best-possible constraints: Ignoring systematic effects 
such as metallicity, dust law, IMF, and choice of SPS 
code, what are the best-possible constraints on the stel- 
lar population parameters of individual galaxies we can 
expe(;t from their SEDs alone? 3) Comparison of SPS 
Codes: How well do the most-used SPS codes describe 
the SEDs of young stellar populations in the rest-frame 
NIR? Including the IRAC data, do any of the codes pro- 
vide better fits to the data? What are the systematic dif- 
ferences in parameters for individual galaxies using vari- 
ous SPS codes, dust laws, and metallicities, and which of 
these effects is most significant? 4) Are passive galaxies 
really passive ? Approximately 50% of the galaxies in our 
sample show no emission lines in their spectra. Do all of 
these galaxies have strongly suppressed star formation, 
or are some dusty star forming galaxies with enough dust 
to obscure the emission-line regions? 

This paper is organized as follows. In § 2 we briefly 
review the observational data used in our analysis. In 
§ 3 we introduce the SPS codes used in our fitting, and 
discuss our fitting and error estimation methods. In § 

4 we evaluate the constraints on stellar population pa- 
rameters when including the rest-frame NIR data. In § 

5 we compare the quality of fits from the various SPS 
codes and compare the systematic changes in parame- 
ters from these codes. We assess the systematic effects 
from the various assumptions of dust law and metallicity 
compared to SPS code in § 6, and in § 7 we examine 
what fraction of the galaxies have strongly suppressed 
star formation based on their SEDs. We conclude with 
a summary in § 8 

Throughout this paper we assume a Clm = 0.3, Cl\ = 
0.7, Ho = 70 km s~^ Mpc"^ cosmology. 

2. DATA 
2.1. The Galaxy Sample 



3 



The 34 galaxies^ used in this study are those with 
NIR spectroscopic observations performed by Kriek et 
al. (2006; 2007), K08. Galaxies were selected for the 

sample based on their K-band magnitude (K < 19.7, 
Vega magnitude) and that they were likely to lie in the 
2.0 < z < 2.7 redshift range based on their photometric 
redshift probability distribution as derived from broad- 
band UBVRIz'JHK photometry. As discussed in K08, 
Mann- Whitney and Kolmorov-Smirnov tests show that 
this sample is representative of the distribution of the 
-2phot; J-K, R-K, and (U-V)rest colors of a mass-limited 
sample at 2 < z < 3. It may be less representative of 
a K-bright subsample of the 2 < z < 3 population be- 
cause the overall K-bright population tends to have a 
lower median redshift and narrower redshift distribution 
than the spectroscopic sample (K08). Small biases like 
these should not affect the analysis in this paper. Ap- 
proximately 20% (7/34) of the sample have spectroscopic 
redshifts 1.5 < z < 2.0. These galaxies may not neces- 
sarily be representative of the stellar populations of a 
mass-limited sample in this redshift range; however, we 
include them in this analysis. 

2.2. Photometric Data 

The UBVRIz'JHK photometry for galaxies in the fields 
SDSS1030, CW1255, HDFSl, and HDFS2, (hereafter, 

the MUSYC fields) is derived from the catalogues pre- 
sented in Gawiser et al. (2006) and Quadri et al. (2007). 
Photometry in the same bands for galaxies in the ECDFS 
field are derived from the catalogues presented in Taylor 
et al. (2008). The IRAC photometry for galaxies in the 
MUSYC fields is presented in Marchesini et al. (2008) 
and we refer to that paper for a detailed discussion of 
the data reduction and catalogue creation. The IRAC 
data for the ECDFS field was obtained as part of the 
SIMPLE survey and is discussed in Damen et al. (2008). 
We have added a systematic error of 10% in quadrature 
to the IRAC photometric errors to account for zero point 
uncertainties and the color-dependent flat fielding errors. 

A few of the galaxies in the K08 sample were located 
close to the edge, or completely out of the region covered 
by the IRAC data. SDSS1030-101 is missing ^.Q^iva and 
5.8/um data, whereas SDSS1030-1839 is missing 4.5/xm 
and 8.0/im data. These galaxies still have data in the 
complementary IRAC bands so we retain them in the 
sample. The galaxies SDSS1030-301 and SDSS1030-1813 
were near the edge of the field and due to the dither pat- 
tern have only a few frames of data in the IRAC bands. 
Although they arc detected, the lack of data make back- 
ground subtraction and cosmic ray rejection difficult so 
we remove them from the sample. Removal of these 
galaxies reduces the sample to 34 of the 36 galaxies pre- 
sented in K08. 

2.3. NIR Spectroscopic Data 

The reduction of the NIR spectroscopic data used in 
our SED modeling was discussed in detail in Kriek et al. 
(2006). Briefly, the data were obtained with the GNIRS 
instrument on Gemini-South using the 0.675"-wide slit 
and the 32 1/mm grating in cross-dispersed mode which 
provided a resolving power of R 1000. Spectroscopic 

^ The original sample contains 36 galaxies; however, two were 
not observed as part of the IRAC program, see § 2.2 



redshifts were determined for 17/34 galaxies that had de- 
tected emission lines. For the remaining 17 galaxies that 
did not have detectable emission lines K08 determined 
a spectroscopic continuum redshift using SED modeling 
of the NIR spectra and broadband photometry. In this 
analysis we use the low-resolution binned spectra con- 
structed by K08. Each data point in the binned spectra 
contains 80 good pixels from the observed spectrum, cor- 
responding to ^ 400A observed frame. Absolute flux cal- 
ibration of the spectra is performed using the broadband 
JHK magnitudes. 

3. SED FITTING AND STELLAR POPULATION 
PARAMETERS 

3.1. Stellar Population Models 

At present, there are several well-tested SPS codes 
pubhcally available. Codes such as PEGASE (Fioc 
& Rocca-Volmerage 1997), Bruzual & Chariot (2003, 
hereafter BC03), Maraston (2005, hereafter M05), Star- 
burstQQ (Vazquez & Leitherer 2005), and Chariot & 
Bruzual (2008, in preparation, hereafter CB08) all pro- 
vide synthetic spectra for SED modeling. 

As pointed out by M05, the PEGASE, Starburst99, 
and BC03 models use similar stellar spectral libraries, 
isochrones, and treatment of the TP-AGB phase of stel- 
lar evolution. Therefore, in our fitting we use the BC03 
models assuming they are representative of this class. 
We also perform fits using the M05 code which uses the 
fuel consumption approach as an alternative to isochrone 
synthesis for modeling stellar evolution off the main se- 
quence. Lastly, we perform fits with the CB08 code 
which is similar to the BC03 code except that it also 
includes updated evolutionary tracks for TP-AGB stars 
from Marigo et al. (2008). Taken together these three 
codes span the different treatments of TP-AGB stars cur- 
rently available. 

3.2. Fitting Method 

The fitting procedure used in our analysis is analo- 
gous to the one outlined in Kriek et al. (2006) and K08. 
The photometric and spectroscopic data are fit to mod- 
els with 34 exponentially declining SFHs (parameterized 
by r) ranging from 0.01 to 20.0 Gyr. We fit 45 different 
ages [t) ranging from 0.01 to 13 Gyr, but only allow those 
that are less than the age of the universe at the redshift 
of the galaxy. For the majority of comparisons in this 
paper we assiimc solar metallicity and a Calzetti et al. 
(2000) dust law^ as our "control" model. The V-band 
attenuation (A„) is allowed to range between and 4 
mag in increments of 0.1. When fitting galaxies without 
emission lines we also adopt z as a, free parameter and 
fit in increments of 5z = 0.01. 

To maintain continuity with the fits presented in K08 
we adopt a Salpeter (1955) IMF. Although it is widely 
accepted that the Kroupa (2001) and Charbrier (2003) 
IMFs are likely to be better choices for a universal IMF; 
we note that the galaxies in our sample are young and 
the main sequence turnoff mass is much larger than the 
mass at which these IMFs differ, therefore, the choice be- 
tween these IMFs does not significantly affect the shape 
of the synthetic SEDs, only the normalization of the 

® A comparison of fits using other dust laws and metallicities is 
presented in § 6 
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Mstar- Given an average age of 1 Gyr, the Mgtar com- 
puted for our sample can be converted to a Charbrier or 
Kroupa IMF simply by multiplying by scaling factors of 
0.57 and 0.63, respectively. Our values of Mgtar are also 
corrected to account for gas recycled to the interstellar 
medium from supernovae and planetary nebula. There- 
fore, our Mstar corresponds to living stars plus stellar 
remnants, but do not include recycled gas. 

For each SPS SED, we find a normalization (X) by 
minimizing the function 

2 1 ' fmodel{T,t,A„,z),i fobs,i) 

where Ndof is the number of degrees of freedom. lobs,i 
is the observed flux in the i*'* filter, 5iobs,i is the error 
in the observed flux, and imodei(T.t.A^,z),i is the flux of 
model SED which is a function of t, i, A^,, and z. We 
then construct a x^-surface as a function of the synthetic 
SPS model parameters t, t, A^,, and z. The best fit 
SPS model is determined by finding the location of the 
minimum on the Xr surface. 

The la confidence intervals in the model parameters 
are determined using the Monte Carlo method suggested 
by Papovich et al. (2001) and K08. In these simulations 
the photometric and spectroscopic data are perturbed 
randomly within their uncertainties and then the SPS 
models are fit again. From 200 simulations per galaxy we 
find the value of that encompasses the Xr from 68% of 
the Monte Carlo simulations. Returning to the original 
X^ fitting surface, locations with x^ below this value are 
considered to be allowable within the la errors. The 
strength of this method is that by using the original Xr" 
surface the correlations between errors in the parameters 
are preserved. 

4. IMPROVEMENT IN STELLAR POPULATION 
PARAMETERS FROM REST-FRAME NIR DATA 

In this section we examine the importance of IRAC 
data for constraining the stellar population parameters 
of 2; ~ 2.3 galaxies from their SEDs without considering 
systematic effects such as the choice of metallicity, dust 
law, IMF and SPS code. We test for potential systematic 
differences in parameters that result from fitting with 
and without constraints on the rest-frame NIR SED, as 
well as how much (if at all) the rest-frame NIR improves 
the uncertainties in the stellar population parameters. 
Throughout this section we use the BC03 r-models with 
solar metallicity, the Calzetti et al. (2000) dust law, and 
a Salpeter IMF as our "control" model for the compar- 
isons. 

4.1. Comparison of SED Fits With and Without IRAC 

Data 

In Table 1 we list the best fit SED parameters and 
la uncertainties determined from fits to the UBVRIz' 
broadband data and NIR spectroscopy (i.e., with- 
out the IRAC data, hereafter we refer to this as 
U^z'-|-NIRspec), as well as those from fits with the 
UBVRIz', NIR spectroscopy and the IRAC data (here- 
after U— >8/um-|-NIRspec). Table 1 also contains the pa- 
rameter (t)sFR, defined as the SFR- weighted mean age of 
the stellar population (see Forster Schreiber et al. 2004). 



The parameter t in Table 1 is the time since the onset of 
SF; however, in r-models galaxies are continually form- 
ing stars and therefore (t)sFR defines the mean age of the 
stellar population that dominates the light of the galaxy. 
Throughout this paper we use (t)sFR as the metric of the 
age of the galaxies. 

K08 also determined SED parameters for this sample 
of galaxies using the U^z'-|-NIRspec data (see their Ta- 
ble 2). Some of our U^z'+NIRspec fit parameters are 
identical to those determined by K08; however, most of 
the parameters and uncertainties have changed by small 
amounts. These minor differences occur for several rea- 
sons. Firstly, our grid of t and t is larger and the spac- 
ing is somewhat different than the K08 grid. Secondly, 
some of the broadband fluxes have changed slightly in our 
newer photometric catalogues. Thirdly, we use a new fit- 
ting code which finds the SED normalizations using a dif- 
ferent algorithm than the K08 code. This leads to small 
numerical differences in the Xr'^- Frequently the differ- 
ence between the minimum Xr of the K08 best fit and 
our best fit is < 0.1%; however, the best fit parameters 
can be significantly different. There are a few galaxies 
that have large differences in parameters; however, the 
la confidence intervals from our fits and the K08 fits still 
span the same range of values for these galaxies. Due to 
these changes we use the parametcirs and uncertainties 
for the U^z'-hNIRspec data determined using our fit- 
ting code for consistency of comparison. 

In Figures 1 and 2 we plot F^ as a function of observed 
wavelength for the photometric data as well as the best 
fitting SEDs. Broadband photometry is plotted in red 
and the binned NIR spectroscopic data is plotted in cyan. 
The SEDs fit with the U^z'-|-NIRspec data are plotted 
in grey, and the SEDs fit with the U^8/xm+NIRspec 
data are plotted in black. The fit parameters are listed 
in the panels of Figures 1 and 2 for ease of comparison. 
The broadband JHK fiuxes have been plotted for refer- 
ence, but were not used in the fitting. 

Examination of Figures 1 and 2 shows that for ^ 50% 
of the galaxies the SEDs fit without the IRAC photom- 
etry are still consistent with the IRAC photometry at 
< la. Comparing the best fit SED parameters we find 
that for 27/34 galaxies (79%) the z, r, t, Ay, M^tar, SFR, 
and (t)sFR determined without IRAC data all agree with 
those fit with IRAC data within the 1 a confidence in- 
tervals. This demonstrates that for most galaxies any 
systematic changes in the fits that occur due to the in- 
clusion of the IRAC data are smaller than the formal 
confidence intervals (we explore this in more detail in § 
4.3). For the remaining 7/34 galaxies (1030-1531, 1030- 
2728, 1256-519, 1256-1967, HDFS2-509, ECDFS-6842, 
ECDFS-12514) there is a mixed range of agreement be- 
tween the black and grey SEDs and the best fit parame- 
ters. 

The largest difference in best fit parameters as well as 
SED shape is for the galaxy 1256-519, which provides 
a qualitative illustration of the importance of the rest- 
frame NIR photometry for constraining stellar popula- 
tion parameters. Without the IRAC data galaxy 1256- 
519 is best fit as a young, dusty star forming galaxy {t = 
0.3 Gyr, A„ = 3.2, SFR ~ 400 Mg yr'^), but including 
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the IRAC data it is best fit as an old'', moderately-dusty 
quiescent galaxy {t = 2.5 Gyr, = 0.9, SFR ~ 3 Mq 
yr^^). As Figure 1 illustrates, there is little difference be- 
tween these two SEDs in the rest-frame optical; however, 
their fluxes in the rest- frame NIR are quite different. A 
young- and- dusty population is much brighter in the rest- 
frame NIR than an old-and-quiescent population. Once 
the IRAC data is included it is clear that the galaxy has 
an SED that is consistent with an old-and-quiescent pop- 
ulation. 

The importance of the rest-frame NIR for distinguish- 
ing young-and-dusty populations from old-and-quiescent 
populations has already been discussed by previous au- 
thors, e.g., Labbe ct al. (2005), Papovich ct al. (2006), 
Wuyts et al. (2007) and WiUiams et al. (2008). 
These studies have suggested that simple color-color di- 
agrams could be an efBcient method for separating these 
types. Galaxy 1256-519 clearly demonstrates how criti- 
cal the rcst-framc NIR is; even with spectroscopic red- 
shifts and high-resolution spectrophotometry near the 
Balmer/4000A break, it can still be difficult to robustly 
distinguish young-and-dusty from old-and-quiescent pop- 
ulations without additional data on the rest-frame NIR 
SED. 

4.2. Improvement in Photometric Redshifts with IRAC 

data 

Accurate determination of the stellar population pa- 
rameters of z ^ 2 galaxies with broadband data requires 
high quality photometric redshifts (zphot)- Using this 
spectroscopic sample K08 showed that the 2;phot's deter- 
mined using the Rudnick et al. (2001; 2003) code, which 
employs the Coleman et al. (1980) and Kinney et al. 
(1996) templates, were systematically overestimated by 
Az/{1 -1-2;)= 0.08, where Az = (zphot - ^spec), and had 
a scatter of 0.13 in the same units. Although random 
errors in ^^phot can be overcome with larger samples, sys- 
tematic errors can pose significant problems, particularly 
when studying the evolution of the stellar mass density 
or luminosity density (see, e.g., K08, van Dokkum et al. 
2006). It has been suggested by some authors that the 
1.6/im bump feature present in the SEDs of evolved stel- 
lar populations could be a useful Zphot indicator (e.g., 
Simpson & Eiscnhardt 1999; Sawicki 2002). This feature 
falls in the IRAC bands for galaxies at 1.5 < z < 3.0 and 
given the large scatter and systematic offset between the 
Zphot and the spectroscopic redshifts (zspec) seen by K08, 
it is worth investigating whether deep IRAC data might 
improve the Zphot of z ^ 2.3 galaxies. We test this by 
determining Zphot using the broadband photometry with 
and without the IRAC data and comparing these to the 

•^■spec* 

The Zphot's are computed using the EAZY photometric 

redshift code (Brammer ct al. 2008). We use the stan- 
dard EAZY vl.O template set which is determined using 
nonnegativc matrix factorization of SED models from 
the PEGASE code (Fioc & Rocca-Volmerange 1997). 
For comparison, the Zphot are fit using three datasets. 

Throughout this paper wc use the term "old" when referring 
to stellar populations that arc > 1 Gyr old. In the local universe 
this age would be considered young; however, at 2 ~ 2.3 it is nearly 
a maximally old stellar population. Therefore, all galaxies in our 
sample are implicitly "young" , and we use the terms "young" and 
"old" in their relative sense. 



the U— >K photometry, the U— >4.5/um photometry and 
the U^8.0/um photometry. The U^4.5/um photometry 
is tested separately from the U^8.0/^m photometry be- 
cause in practice the 5.8/im and 8.0/im bands arc often 
avoided when fitting Zphot for large samples of galaxies 
(e.g., Brodwin et al. 2006; Marchesini et al. 2008) be- 
cause at present most template sets do not include the 
PAH features that fall in these bandpasses at z < 0.7. 
Wc omit the galaxies 1030-101 and 1030-1839 when com- 
paring the Zphot because they only have photometry in 
two of the IRAC bands. 

In Figure 3 we plot Zphot vs. z^pec for the three 
datasets. For galaxies without emission lines the Zgpcc is 
the best-fit redshift using the U^8/im+NIRspcc data in 
§ 4.1, including the appropriate error bar. Following K08, 
galaxies in the ECDFS field, which has much lower S/N 
photometry in the JHK bands than the other MUSYC 
fields are plotted with open grey circles. Comparing the 
•s^phot fit with the U^K photometry to the Zgpcc we find 
a systematic offset of 0.17 in Az/ (1 -I- z) and a scatter of 
0.12. This offset and scatter are similar, but not identi- 
cal to the ones measured by K08; however, they used the 
Rudnick et al. (2001; 2003) code rather than EAZY, so 
we do not expect perfect agreement. Most of the system- 
atic offset in the Zphot is caused by the ECDFS galaxies 
(hereafter the "wide" sample, following the convention 
from K08). Removing those from the sample (hereafter 
the "deep" sample) the systematic offset is 0.00 and the 
scatter is 0.05, comparable to the offset and scatter of 
0.03 and 0.08 measured by K08. Wc list the offsets and 
scatters measured for all datasets in Table 2. 

When the Zphot are computed with the U^4.5/Lim pho- 
tometry, there is no improvement in the offset and scat- 
ter. Surprisingly, adding the 3.6/im and 4.5/im data does 
not even improve the Zphot of the wide sample, which has 
the lowest S/N ncar-IR photometry. When we fit with 
the entire U— »8/um data the offset and scatter in the over- 
all sample improve to 0.14 and 0.11, respectively, and this 
overall improvement is caused primarily by improvement 
in the wide sample where the systematic offset is reduced 
by ^0.1 in 5z. For the deep sample the offset and scatter 
are 0.02 and 0.05, comparable to what is found with just 
the U— >K photometry. 

These comparisons show that the IRAC photometry 
docs modestly improve the Zphot of galaxies at z ^ 2, but 
it does so only for galaxies with marginal S/N photom- 
etry in the optical and NIR bandpasses, and only when 
all four IRAC channels are used. For galaxies with good 
S/N in the optical and NIR bandpasses, there is no evi- 
dence for an improvement in the photometeric redshifts 
when IRAC data is included. This is best illustrated 
by comparing the Zphot of the deep sample, fit without 
the IRAC data to the wide sample, fit including all four 
IRAC channels. For the deep sample without IRAC data 
the offset and scatter are 0.00 and 0.05, respectively, and 
for the wide sample with IRAC data they are 0.35 and 
0.17, respectively. Clearly deep NIR data arc more use- 
ful than IRAC data when determining the Zphot's of z 
2.3 galaxies. This suggests that for this sample, the con- 
straints on the Zphot are dominated by the location of the 
Lyman, Balmer or 4000A breaks in the SEDs, and that 
constraints from the 1.6/um bump are much weaker. The 
Lyman, Balmer and 4000A breaks are sharper features 
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Fig. 1. — SEDS fit using BC03 models without IRAC data (grey SED) and with IRAC data (black SED). The red points are the 
broadband photometry and the cyan points are the binned NIR spectroscopy. 



than the 1.6/xm bump, and the optical and NIR bands 
are closer spaced in wavelength than the IRAC bands, 
and therefore it is not surprising that this is the case. 

It is worth pointing out that although the ^phot of the 
1.6 < z < 2.9 spectroscopic sample are not significantly 
improved with the IRAC data, this does not necessar- 
ily mean that IRAC data will not improve Zphot for a 
different sample. Our K-selected sample is comprised 
primarily of massive galaxies, which are frequently red 



with strong Balmer/4000A breaks. Blue galaxies with a 
weak Balmer break and the Lyman break blueward of the 
observed U-band (z < 2.5) may see larger improvements 
in 2phot when IRAC data is included. 

4.3. Reduction of Systematic Errors in Stellar 
Population Parameters with IRAC Data 

Examination of the SEDs in § 4.1 showed that 
when the IRAC data were included in the fits to the 
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Fig. 2. — As Figure 1. 

U^z'+NIRspec data, the stellar population parameters 
of most galaxies remain consistent within the la uncer- 
tainties. Although consistent, the differences may be sys- 
tematic, which could change the mean stellar population 
parameters of the sample. In this section we quantita- 
tively check for systematic differences in the mean pa- 
rameters with the inclusion of the IRAC data. 

In Table 3 we list the mean value of the parame- 
ters determined with different subsets of the data com- 



pared to the mean value determined using the full 
U^8/im+NIRspec data set. Specifically, we compare pa- 
rameters from fits to the broadband photometry and NIR 
spectroscopy without the IRAC data (U^z'-|-NIRspec, 
i.e., the same data as used in the fits by K08), broad- 
band photometry including the IRAC data (U^8/Ltm), 
and broadband photometry without IRAC data (U^K). 
When fitting the broadband data we leave redshift as a 
free parameter in order to preserve systematic differences 
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Fig. 3. — Top panel: EAZY ^p^ot calculated using the U—>K pho- 
tometry vs. Zspec- Galaxies in the ECDFS, which has much poorer 
S/N JHK photometry than the MUSYC deep fields are plotted as 
open grey circles. Middle panel: Same as top panel but for Zphot '^ 
determined with U— >4.5//m photometry. Bottom Panel; Same as 
top panel but for ZpjjQj's determined with U— >8.0/im photometry. 
The IRAC data docs not improve the Zphot °f galaxies with good 
S/N JHK photometry. It does reduce the systematic overestimate 
of Zphot galaxies with poorer S/N JHK photometry, but only 
by 0.1 in 5z, and only when all four IRAC channels are used. This 
shows that high S/N JHK photometry is significantly more valu- 
able than IRAC data when determining the Zphot '^^ galaxies at 
z ~ 2.3. 

that may result from either random or systematic errors 
in Zphot- The uncertainties listed in Table 3 are stan- 
dard errors of the mean {a/\fN). More details of the 
uncertainties and histograms that show the shape of the 
distributions are presented as an Appendix. 

Table 3 shows that changes in the mean values of Mgtar , 
"T, (^)sFRj At,, and SFR determined without the IRAC 
data (U^z'-|-NIRspec) compared to those determined 
with the IRAC data (U— >8/xm-|-NIRspec) are all consis- 
tent within the statistical accuracy achievable with our 



sample, approximately 10%. This demonstrates that not 
only are the parameters determined by K08 consistent 
with those determined with the IRAC data, but that 
none of the differences serve to change the mean param- 
eters of the sample in a systematic way. 

The mean values of most parameters determined us- 
ing the broadband data are also consistent with those 
determined using the U-^-S^m-l-NIRspec data. The only 
exception is Mstar, which appears to be overestimated 
by a factor of ~ 1.3 (0.11 dex) when only U^K pho- 
tometry is used in the fitting. This systematic overes- 
timate was also seen by K08 when they compared the 
Mstar from fits to the U— >K data to those from the 
U^z'+NIRspec data. Given that Mgtar determined with 
both the U^8/xm-|-NIRspec and U^z'-|-NIRspec data 
are consistent, our fits support the K08 conclusions. As 
discussed in K08, the systematic effect on Mgtar can lead 
to an overestimate of the number of massive galaxies at 
high redshift, and hence an underestimate of the evolu- 
tion of the stellar mass density. 

Interestingly, it appears that the systematic overesti- 
mate of Mstar is removed when the IRAC data is in- 
cluded in the SED fits to the broadband data. Eisner 
et al. (2008) also saw a systematic reduction in Mgtar 
when including IRAC data in SED fitting of broadband 
photometry of z ~ 2 galaxies. Their Mstar was system- 
atically reduced by a factor of 1.6, which is similar to 
our factor of 1.3. Shapley et al. (2005) and Wuyts et 
al. (2007) made the same comparison; however, neither 
found that adding IRAC data caused a significant sys- 
tematic change in the mean Mstar of their sample. Shap- 
ley et et al. (2005) did suggest that the Mgtar of a sub- 
sample of their galaxies was systematically reduced when 
IRAC data was included in the fitting. 

Why does including the IRAC data remove the sys- 
tematic overestimate of Mgtar in owi sample? In § 4.2 we 
showed that the z^hot were systematically overestimated 
for galaxies in the ECDFS field, and that the IRAC data 
reduced this systematic by ~ 0.1 in Az/(1 -I- z). If we 
remove the ECDFS galaxies and compare the mean val- 
ues of Mstar from the U^K photometry again, we find 
that the systematic overestimate of Mgtar is reduced to a 
factor of 1.08 ± 0.10. This suggests that it is primarily 
the systematic overestimate of the Zphot, caused by low 
S/N JHK photometry that is responsible for the overes- 
timate of Mstar- Given that Wuyts et al. (2007) already 
had high S/N JHK photometry, this may explain why 
they saw no change in Mgtar when the IRAC data was 
included in the fitting. Both the Eisner et al. (2008) and 
Shapley et al. (2005) studies use the same redshifts for 
their fits with and without IRAC data, which suggests 
that their differences may be more data-specific, i.e., de- 
pend on the combination of broadband filters. Shapley 
et al. (2005) suggest it may be caused by contamination 
of their observed K^-band fluxes from Ha emission. 

In summary, our comparison shows that when compar- 
ing stellar population parameters determined with var- 
ious permutations of the data, the only parameter that 
suffers a systematic bias at the > 10% level is Mgtar 
when determined using the U^K photometry. This bias 
can be removed by including either NIR spectroscopy 
or IRAC data in the fitting because both improve the 
^phot- None of the other parameters, r, (t)sFR, A„, or 
SFR show evidence for systematic differences at > 10% 
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even without using IRAC data or NIR spectroscopy in 
the fitting. 

4.4. Improvement of Uncertainties in Stellar Population 
Parameters with IRAC Data 

Although the IRAC data do reduce the systematic er- 
rors in Zphot and Mgtar for our sample, there are no 
changes in the mean values of the other stellar popu- 
lation parameters. This suggests that the IRAC data 
do not improve the accuracy for the majority of stellar 
population parameters. In this section we examine if in- 
cluding the IRAC photometry improves the uncertainties 
in the parameters determined for individual galaxies. 

A simple test is to compare the parameter uncertain- 
ties computed using the Monte Carlo method. While 
potentially informative, such a comparison requires the 
uncertainties themselves to be robust, even in the case 
of sparse data. A direct way to test for an improvement 
is to measure the distribution of parameters computed 
without IRAC data relative to the distribution computed 
with IRAC data. The rms scatter of this distribution is 
a metric of the average uncertainty in a parameter for 
the sample, independent of the method used to estimate 
the uncertainties for individual galaxies. 

Measuring a relative distribution of parameters re- 
quires a comparison sample. Since we do not have in- 
dependent knowledge of the parameters without mea- 
surement uncertainty, we will assume that the parame- 
ters determined using the U^8/im+NIRspec data are 
likely to be the most precise, and compute the stan- 
dard deviation (a) of parameters determined without the 
IRAC data relative to these. A more detailed discussion 
of these comparisons, as well as histograms of the dis- 
tributions are presented in the Appendix. In Figure 4 
we graphically summarize the result of the comparisons 
for the parameters Mgtar, t, (i)sFR, A^, and SFR, com- 
puted using the U^K, U^8/im, U^z'-|-NIRspec, and 
U— >8/xm-|-NIRspec data. 

The left panel shows the cr of the logarithm of the pa- 
rameter distributions computed using each subset of the 
data plotted as a function of the a of the parameter dis- 
tribution from the U^8/um+NIRspec data. In order to 
compare parameters with different units we use the frac- 
tional scatter in each parameter ((j/mean). Reading from 
left to right along the X-axis of Figure 4 shows which pa- 
rameters are determined with the best precision using the 
U^8/xm+NIRspec data. 

The Y-axis of Figure 4 shows the logarithm of the scat- 
ter in parameters determined using the U^K, U^8/im, 
and U^z'-|-NIRspec data which are plotted as green di- 
amonds, blue squares, and red circles, respectively. The 
one-to-one relation is plotted as a dotted line. If a data 
point lies on that line it indicates that the precision of 
that parameter, with the given subset of data, is as good 
as the precision of that parameter attainable with the 
entire dataset. 

Figure 4 shows some interesting trends, perhaps the 
most unsurprising of which is that the uncertainties in 
parameters determined using the U— >K data are always 
the largest. Adding either the NIR spectroscopy, the 
IRAC data, or the combination of the two significantly 
improves the uncertainties in all parameters. If we com- 
pare the uncertainties in parameters determined with the 
U— >8/um data (blue squares) to those determined with 



the U— >z'-|-NIRspec (red circles) we can assess which pa- 
rameters are most improved with each type of data, as 
well as the size of the improvement. 

Examining the constraints from the V^Sfim data, it 
is clear that the IRAC data is most useful for constrain- 
ing the Mstar, A„, and SFR. Interestingly, adding the 
NIR spectroscopy in combination with the IRAC data 
does not further improve these parameter estimates. The 
U— >z'-|-NIRspec data shows that the NIR spectroscopy 
is most useful for constraining the (t)sFR and t of the 
models; however, the uncertainty in both of these pa- 
rameters can be further improved by including the IRAC 
data. The NIR spectroscopy also improves the Mgtar sig- 
nificantly, but including the IRAC data makes the con- 
straints only ^ 5% better. 

Why do the different types of data affect the param- 
eters in these ways? How can IRAC data improve SFR 
estimates, when it primarily traces the old stellar popula- 
tion? Although young-and-dusty and old-and-quiescent 
systems are difficult to distinguish with data that cov- 
ers only the rest-frame UV to optical SED, they are 
separable once rest-frame NIR data is available (e.g., § 
4.1, Labbe et al. 2005; Williams et al. 2008). This is 
not because the NIR wavelength range itself is partic- 
ularly valuable, but because it completes the coverage 
of the part of a galaxy's SED that is dominated by light 
from stellar photospheres. In turn this improves the con- 
straints on the A^ , because A^ changes the broad shape 
of the stellar SED. In the r- models the SFR is not a 
quantity that is fit for, instead it is inferred from the 
number of e- folding times (t/r), scaled by the Mstar, and 
corrected for the A„. Improved constraints on A^ lead 
to better SFRs, and hence IRAC data actually improves 
estimates of SFRs from SED fitting. 

The NIR spectroscopy adds both a Zgpcc, as well 
as high-resolution information on the SED near the 
Balmer/4000A break. As shown by K08, it is the high 
resolution information near the Balmer/4000A break 
that drives the improvement in both t and r. Inter- 
estingly, the IRAC data still helps to improve the con- 
straints on these parameters further. Most likely this 
is because it traces the old stellar population, which in 
tandem with the optical data constrains the ratio of old 
stars to young stars, and hence the number of e-folding 
times, which is directly related to both t and r. 

Perhaps the most remarkable result from Figure 4 is 
that for our sample, the parameters Mgtar and SFR are as 
well-determined from the U— *8/zm photometry and Zphot, 
as they are from the U^8/im+NIRspec data and Zspec- 
This suggests that deep broadband photometry that ex- 
tends into the rest-frame NIR can potentially provide 
unbiased high-quality estimates of the Mgtar and SFR 
for massive galaxies at 2 < < 3. NIR spectroscopy is 
extremely valuable for getting accurate rest-frame colors 
using Zspec (e.g., K08), identifying AGN (e.g., Kriek et 
al. 2007), and to provide an independent check on SFRs 
from emission line fluxes (e.g., Erb et al. 2006); however, 
it does not significantly improve the constraints on Mgtar 
or SFR from the SED fits if IRAC data is available. 

4.5. How Well Can We Constrain Stellar Populations 
Parameters with SEDs? 

Given that our current sample of galaxies has the best- 
constrained SEDs of 2; 2.3 galaxies currently available, 
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we can estimate the best-possible constraints on stel- 
lar population parameters that can be determined from 
SED fitting ignoring systematic effects such as choice 

of metallicity, IMF, and dust law. Assuming that the 
BC03 models with solar metallicity, a Salpcter IMF, and 
a Calzetti et al. (2000) dust law create an appropriate 
set of template SEDs, the X-axis of Figure 4 shows that 
the SED fitting best-constrains the parameters Mgtar, A„, 
and (t)sFR- Using the U^8/xm+NlRspec data, Mgtar can 
be determined to ± 0.12 dex, (t)sFR to ± 0.26 dex, and 
A„ to ± 0.3 mag. Interestingly, when only broadband 
U^8/im data is used, the constraints on these parame- 
ters are similar, ± 0.11 dex, ± 0.36 dex, and ± 0.3 mag, 
for Mstar, (i)sFR, and A^,, respectively. 

The r and SFR, are not as well-determined as the other 
parameters. With the entire U^8;um+NlRspec data set 
the scatter in r and SFR are factors of 2.5 and 2.8, respec- 
tively. The uncertainty in the SFR with the broadband 
data alone is still a factor of 2.8; however, the uncertainty 
in T is much worse, a factor of 4.6. . 

Overall, this exercise shows that the best constrained 
parameters of individual galaxies in our sample from 
SED fitting are Mgtar, (^)sfr , and A„. As we will demon- 
strate in the § 5, these uncertainties are now small enough 
that they are comparable to the level of the systematic 
errors caused by the uncertainty in SPS code, metallicity, 
dust law and IMF. 

5. COMPARISON OF SPS CODES 

In this section we fit the full U-^8/xm+NlRspec photo- 
metric dataset to a grid of r-models from the SPS codes 
of BC03, M05, and CB08. We examine the overall quality 
of the fits based on the Xr statistic as well as the residuals 
from the mean SEDs to test whether a particular set of 
models provides a better description of the data. Given 
the different treatment of the TP-AGB phase of stellar 
evolution in the M05 and CB08 models compared to the 
BC03 models and the larger rest-frame NIR fluxes this 
produces, we compare fits with and without the IRAC 
data as a test of how well the models describe the SEDs in 
the rest-frame NIR. We also examine the systematic dif- 
ferences in the stellar population parameters determined 
with the different codes and their implications for studies 
of the stellar populations of distant galaxies. Throughout 
this comparison we assume solar metallicity, a Salpeter 
IMF, and the Calzetti et al. (2000) dust law as a "con- 
trol" model. 

5.1. Comparison of Models When Excluding IRAC 

Our first comparison of the models is made using the 
U^z'-l-NlRspec data which effectively spans the UV 
through optical wavelength range for the galaxies in our 
sample. The parameters determined from these flts are 
listed in Table 1 . In the top panels of Figure 5 we plot 
the x^'s of the fits to the M05 and CB08 models against 
those from the fits to the BC03 models. Figure 5 shows 
that there is a large range of Xr values for the fits; how- 
ever, the majority of the systems, 26/34 (76%) are rea- 
sonably well described by the models, having xi < 2. 
This shows that even with the high-resolution photo- 
metric data available from the NIR spectroscopy, simple 
T-models from all three SPS codes can still successfully 
reproduce the SEDs of the majority of 2; ^ 2 galaxies. 

Comparison of the Xr values shows that the BC03 mod- 



els have the lowest Xr for 9/34 galaxies (26%), the M05 
models have the lowest Xr 15/34 galaxies (44%), and 
the CB08 models have the lowest Xr for 10/34 galax- 
ies (29%). This demonstrates that within our sample 
there is no statistically significant evidence that one of 
the models provides the best fits more frequently than 
the others. Summing the xf values for all 34 galaxies 
results in total Xr's of 57.72, 60.54, and 58.02 for the 
BC03, M05, and CB08 models, respectively. These im- 
ply an average xl of 1.70 ± 0.03, 1.78 ± 0.03, and 1.71 
± 0.03 for the BC03, M05, and CB08 models, respec- 
tively, where the quoted errors are standard errors of the 
mean. This suggests that the BC03 and CB08 models 
may provide modestly better fits on average to the rest- 
frame UV through optical SEDs than the M05 models; 
however, the majority of this difference is caused by a 
few galaxies that are fit much better by the BC03 and 
CB08 models (e.g., galaxy HDFS2-509 and 1256-1967) 

It is worth noting that in Figure 5 the Xr values of the 
BC03 and CB08 models are nearly identical for all galax- 
ies, whereas there is a larger scatter in the x^'s between 
the BC03 models and the M05 models. This confirms 
that the SEDs in the BC03 and CB08 models appear to 
be nearly identical in the rest-frame UV through optical 
wavelength range. In § 5.2 we will show that the fits 
to the rest-frame UV to NIR using the BC03 and CB08 
models are quite different; and the current comparison 
makes it clear that those differences must be driven ex- 
clusively by the rest-frame NIR part of the SED. 

Overall, based on the average Xr, the BC03 and CB08 
models perform slightly better than the M05 models at 
fitting the rest-frame UV to optical SEDs, but the dif- 
ference is small (^ 5% in xl) and only significant at 
1.5cr within our sample of 34 galaxies. 

5.2. Comparison of Models When Including IRAC 

We next compare how well the models describe the en- 
tire rest-frame UV to NIR SED of the galaxies by fitting 
the U^8/^m+NIRspec data. The best fit SEDs from 
each SPS code are plotted in units of F,y vs. A in Figures 
6 and 7. The black line, green line, and blue line are the 
best fits using the BC03, M05, and CB08 models, respec- 
tively. The best fit stellar population parameters, their 
uncertainties, and the Xr of the best fits are listed in Ta- 
ble 1. We plot the Xr's of the fits to the M05 and CB08 
models against those from the fits to the BC03 models 
in the bottom panels of Figure 5. 

Comparison of the (Xr)'s in the top and bottom panels 
Figure 5 shows that the variance in Xr's between mod- 
els is significantly larger when the rest-frame NIR is in- 
cluded in the fitting; however, the frequency with which 
each model provides the best fit is still similar. The BC03 
models have the lowest Xr for 11/34 (32%) of the galax- 
ies, the M05 models have the lowest Xr for 12/34 (35%) 
of the galaxies, and the CB08 models have the lowest Xr 
for 11/34 (32%) of the galaxies. 

Totaling the Xr for all fits with each SPS codes gives 
values of 62.80, 62.11, and 64.72 for the BC03, M05, and 
CB08 models, respectively. Converting the total Xr's to 
a iXr) P'2r galaxy for our sample gives values of 1.85 
± 0.03, 1.83 ± 0.03, and 1.90 ± 0.03 for the BC03, 
M05, and CB08 models, respectively, where the errors 
are the standard error of the mean. Interestingly, when 
the rest-frame NIR is included, the M05 models describe 



11 



1.5 



X 
0) 

2, 

0) 

CO 

15 
c 
g 

CO 



^ 0.5- 



_2 

o 
w 

< 



0.0 



—I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r- 



■ U^8/im 

• U^z' + NIRspec 
□ U^8/im+NIRspec 



SFR 
▲ 



Tau 



Av 
▲ 



_l I I I I I I I l_ 



Age 



_l I I I I I I I L_ 



_l I I I I I I I L_ 



_l I I I I I I I L_ 



0.1 



0.2 0.3 0.4 

Absolute Fractional Scatter U-*8^m + NIRspec (dex) 



0.5 



Fig. 4. — The logarithm of the absolute fractional scatter in parameters, Log(a-/mean), determined with various subsets of the data 
(U— >K, U— >8/im, U— >z'+NIRspec) plotted versus the logarithm of the absolute fractional scatter determined with the entire dataset 
(U— »8/im+NIRspec). The dotted line has a slope of unity. Reading along the X-axis, the best constrained parameters are Mgtan ^v, and 
(*}sFRi which have rms scatters of ~ 0.12 - 0.25 dex. The r and SFR are not as well constrained from the SED alone, and have a scatter 
of a factor of ^ 2-3. 



the SEDs slightly better than both the BC03 and CB08 
models on average. The difference is not significant when 
comparing the BC03 and M05 models; however, it is sig- 
nificant at ^ 1.5(7 when comparing the M05 and CB08 
models. Maraston et al. (2006) found similar results us- 
ing a sample of 7 galaxies at 1.4 < z < 2.7 selected from 
the GOODS survey. For their sample the (Xr) was 1.38 
± 0.10 for fits with the M05 models, and 1.51 ± 0.10 
for fits with the BC03 models, where we have computed 
the error bars based on the standard error of the mean. 
Their results also show a small, but not statistically sig- 
nificant, advantage for the M05 models compared to the 
BC03 models when the entire rest-frame UV to NIR SED 
is fit. 

If the M05 models had the largest (xr) when fitting the 
rest-frame UV to optical SEDs of the galaxies, but the 
smallest (Xr) when the rest-frame NIR data is included, 
it suggests that they may describe the rest-frame NIR 
for these galaxies better than both the BC03 and CB08 
models. We test this by plotting the Xr values obtained 
when fitting using only the U^z'-|-NIRspec data versus 
those obtained when fitting the U^z'-|-NIRspec-|-IRAC 
data for all three SPS codes in Figure 8. Galaxies that 
have much poorer fits when the rest-frame NIR is in- 
cluded will lie to the right of the dotted line in Figure 8. 

Before we discuss the Xr values from fits with and with- 
out the rest-frame NIR data, it is worth noting that the 
fraction of massive galaxies that host an AGN at ^ ~ 2 
is much higher than at lower redshift (e.g. Kriek et al. 
2007, Daddi et al. 2007). The presence of a dusty AGN 



can cause significant emission at rest-frame NIR wave- 
lengths (e.g., Donley et al. 2008, and references therein), 
and therefore can "contaminate" the stellar SED in the 
rest-frame NIR. Such systems are now being identified at 
high redshift using the presence of excess emission at A > 
1.6/im, rest- frame. Galaxies with a dusty AGN are fre- 
quently referred to as power-law galaxies (PLGs) because 
the strong NIR emission at A > 1.6/im from the AGN 
overpowers the stellar emission removing the turnover 
in the SED at 1.6/im from the peak of the stellar emis- 
sion, and making the SED appear more like a power-law. 
Given that the SPS models do not include an AGN com- 
ponent, we expect any galaxies identified as PLGs will 
be fit much more poorly when the rest-frame NIR is in- 
cluded. Poor fits to these galaxies in the rest-frame NIR 
probably does not reflect a deficiency with the models 
and therefore we need to identify any PLGs within our 
sample. 

We identified PLGs by fitting the IRAC data for each 
galaxy to a power-law of the form (x v°' . Galaxies fit 
with a > -0.5 are considered PLGs, and therefore likely 
candidates for hosting a dusty AGN (e.g., Alonso Her- 
rero et al. 2006, Donley et al. 2007). The properties of 
these galaxies and their overlap with emission-line AGN 
(Kriek et al. 2007) will be presented in more detail, and 
including MIPS 24/ini observations, in a future paper. 
Here we only discuss them as they are relevant to SED 
fitting in the rest-frame NIR. 

The PLGs are plotted as open red circles in Figure 
8. Interestingly, not only do most of the PLGs have 
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Fig. 5. — Top Panels: Comparison of Xr values for fits without the IRAC data included using models from the M05 code vs. the BC03 
code (left) and the CB08 code vs. the BC03 code (right). The (Xr) values arc 1.70 ± 0.03, 1.78 ± 0.03, and 1.71 ± 0.03 for the BC03, 
M05, and CB08 models, respectively. This suggests that on average the BC03 and CB08 models may provide slightly better fits to the 
rest-frame UV to optical of young galaxies than the M05 models. Bottom Panels: Comparison of Xr values for fits including the IRAC 
data using models from the M05 code vs. the BC03 code (left) and the CB08 code vs. the BC03 code (right). The (xj) values are 1.79 ± 
0.03, 1.77 ± 0.03, and 1.85 it 0.03 for the BC03, M05, and CB08 models, respectively. Based on this compajrison, all three codes perform 
equally well at fitting the rest-frame UV to NIR SEDs of young galaxies. 



significantly poorer fits in the rest-frame NIR than the 
average galaxy, they also have signifieantly poorer fits 
in the rest-frame UV to optical than the average galaxy. 
This suggests that if these systems do have a dusty AGN, 
the AGN may also emit radiation blueward of rest-frame 
1.6/v,m and "contaminate" this part of the SED as well. 

Excluding the PLGs, we can see that for the M05 mod- 
els, very few galaxies have significantly higher Xr values 
when the rest-frame NIR data are included in the fit. 
For the BG03 and CB08 models, there are some non- 
PLG systems where the fits are much worse when the 
rest-frame NIR data is included; however, there are only 
a handful of such galaxies. The most discrepant galaxies 
in the BC03 models are 1030-1531, 1256-519, and 1256- 
1967. Galaxy 1256-519 has a poor fit with the IRAC 
data in all models; however, it is notable that 1030-1531 
and 1256-1967 have more flux in the rest-frame NIR than 
predicted by the BC03 models and are better described 
by the M05 and CB08 models. Interestingly, these two 
galaxies also have similar ages (0.47^q'q7, and 0.45^q''|'5 
Gyr) from the M05 fits, and are right in the middle of 
the age range were emission from TP-AGB stars is near 
the maximum for a range of metallicities (e.g. M05) sug- 
gesting that the different treatment of these stars in the 
M05 and CB08 models is the cause of the improved fit. 
Still, it is possible that these two galaxies could be mod- 
eled using a composite of young and old bursts with the 
BC03 code (e.g., Yan et al. 2004), and with only two 



non-PLG galaxies in the sample that would require such 
modeling to explain the SED, we suggest that there is 
only mild evidence at best from the Xr's that the SEDs 
from the BC03 code require significant changes in the 
rest-frame NIR in order to match the observed SEDs of 
our galaxies. 

This can also be seen in Figure 9 where we plot the 
residuals from the SEDs fit using models from the three 
SPS codes, corrected to the rest-frame. The solid line 
represents a running average of the 30 nearest points. 
The running average suggests that the BC03 models may 
underpredict the mean flux of the average galaxy in the 
wavelength range 1.0 - 1.8/um; however, the difference is 
at most, 5-10%, conceivably within the allowable range 
of the uncertanties in the IRAG photometry. The resid- 
uals from the M05 and CB08 models do not show the 
same trend, suggesting that they may describe the rest- 
frame NIR SEDs of these galaxies slightly better than 
the BG03 models. 

In summary, these comparisons show a few key results. 
Most importantly, based on the frequency of lowest Xr 
and the (x^) for all galaxies, we conclude that there is 
no significant evidence that any of the three SPS codes 
we tested describe the complete rest-frame UV to NIR 
SEDs of massive z ^ 2.3 galaxies significantly better than 
the others. Comparison of the (Xr)'s from fits with and 
without the rest-frame NIR data suggest that the BC03 
and CB08 models may describe the rest-frame UV to op- 
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tical SEDs of the galaxies slightly better than the M05 
models, but that the latter may fit the rest-frame NIR 
better than the former. Still, the differences in the (xr) 
between the models arc always < 5% and arc only sig- 
nificant by ~ 1.5ct at most. The residuals from the fits 
also seem to show that the M05 and CB08 models de- 
scribe the rest-frame NIR slightly better than the BC03 
models, but again, the difference is at most 5-10%, and 
within the photometric errors. These comparisons sug- 
gest that it will be extremely difficiilt to refine the models 
based on the observed SEDs of yoimg stellar populations 
at high redshift. Mostly likely such constraints will have 
to come from other avenues such as NIR spectroscopy 
of young stellar populations in the local universe (e.g., 
Riffcl ct al. 2008) or improved theoretical understanding 
of post main sequence stellar evolution (e.g., Marigo et 
al. 2008), or both. 

5.3. How do best-fit parameters vary between SPS 
codes? 

Although the three codes describe the overall shape 
of the SEDs of the galaxies equally well, the parameters 
of the best fit SEDs from the codes imply significantly 
different stellar populations for the same galaxies. In 
Figure 10 we plot the SED parameters determined from 
the U^8/xm+NIRspcc data using the M05 and CB08 
models against those determined using the BC03 models 
as solid circles and open diamonds, respectively. 

5.3.1. Differences in Stellar Populations Between MOB and 
BC03 Models 

The mean and median ratio of stellar population pa- 
rameters determined with the M05 and BC03 models are 

listed in Table 4. The most significant difference between 
these models is in the median Mgtar- From Figure 10 we 
can see that the median Mgtar from the M05 models is 
0.63 that of the BC03 models, and appears to be a sys- 
tematic offset, independent of the mass of the galaxy. 
If we divide the galaxies into two samples, those with 
detected emission lines in the NIR spectrum (hereafter 
EL-galaxies) and those without detectable emission lines 
(hereafter NEL-galaxies) the systematic change in Mgtar 
for both of these groups is 0.63 and 0.69, respectively, 
nearly identical to the median difference of 0.63 for the 
whole sample. This offset in Mgtar is similar to those 
measured by Maraston et al. (2006) and Wuyts et al. 
(2007) who saw ratios of 0.58 and 0.72, respectively. 

Comparing the best fit values of r between the models 
shows that the M05 models prefer slightly lower values, 
they have a median r that is 0.75 that of the BC03 r. 
This difference is driven mostly by the EL-galaxies which 
have a T of 0.66 times the BC03 value; whereas the NEL- 
galaxies have a r 0.83 times the BC03 value. This sug- 
gests that in the M05 models, the SFH of EL-galaxies 
is burstier than for the BC03 models; however, we note 
that r is one of the poorest constrained parameters when 
fitting the SEDs alone. 

If we compare the (t)sFR of the M05 models to the 
BC03 models we find that they are lower by a factor of 
0.65. Again, this ratio is similar to the ratios of 0.58 
and 0.51 found by Maraston et al. (2006) and Wuyts et 
al. (2007), respectively. However, unlike Mgtar, the age 
differences are particularly pronounced between the EL- 
and NEL-galaxies. The ratio of ages of the NEL-galaxies 



between the M05 models and BC03 models is 1.03, but 
for the EL-galaxies the ratio is 0.42. The age of the stel- 
lar population is closely tied to the SFH, as the shape 
of the SED is primarily determined by the number of e- 
folding times of the SFH, which is effectively the ratio of 
t/r. Given that there are differences in r between the 
models, a similar trend in t and (f)sFR is expected in 
order to preserve the number of e-folding times. 

Due to these differences in t and (t) sfr we also expect 
a significant difference in the SFRs for EL- and NEL- 
galaxies. The median ratio of SFRs is 1.12; however, it 
is 1.58 and 0.33 for EL- and NEL-galaxies, respectively 
(the clear trend can also be seen in Figure 11). In the 
M05 models NEL-galaxies have significantly less SF than 
the BC03 models would predict, whereas EL-galaxies are 
more active than the BC03 models would predict, lead- 
ing to a more bimodal galaxy distribution, if galaxies are 
classified by SFR. 

The final parameter we consider is the A^, of the galax- 
ies. The median AA^ for the sample (where AA^ = 
Aw,M05-A„,_BC03) is 0.0; however, AA^, = 0.1 for the EL- 
galaxies, and AA„ = -0.4 for the NEL-galaxies. The 
M05 models require far less dust to fit the NEL-galaxies' 
SEDs, and to first order this explains why the NEL- 
galaxies have smaller Mgtar when fit with the M05 mod- 
els, despite the fact that they have similar r and (t)sFR • 
Maraston et al. (2006) also showed that when fitting 
without dust the M05 models produce fits with signifi- 
cantly lower Xr than the BC03 models, suggesting they 
require less dust to reproduce the SEDs of their sample 
as well. 

In summary, galaxies fit with the M05 models have a 
median values of 0.63-Mstar, 0.65-(f)sFR , and 0.75-r with 

similar AyS and SFRs when compared to galaxies fit with 
the BC03 models, consistent with earlier studies. How- 
ever, breaking the sample up into EL- and NEL-galaxies 
shows more interesting trends. In the M05 models EL- 
galaxies have shorter timescales for SF and younger ages 
and this leads to SFRs that are higher by a factor of 
1.58. Conversely, for NEL-galaxies the values of t and 
(i)sFR are fairly similar between the M05 models and 
the BC03 models; however the A^ is 0.4 mag less in 
the M05 models, and the SFRs are only 0.33 of the 
BC03 value. Overall, the M05 models produce a more 
strongly bi-modal population of galaxies at z ^ 2.3. EL- 
galaxies have a shorter-timescale for SF, and larger SFRs; 
whereas NEL-galaxies have less dust and are more qui- 
escent. 

5.3.2. Differences in Stellar Populations from the CB08 

Models 

Parameters determined from the CB08 models also 
have marked differences from those determined with the 
BC03 models. Similar to the M05 models, Mgtar from 
the CB08 models is lower by a factor of 0.74. Dividing 
these into the EL- and NEL-galaxy categories we find 
that for the EL-galaxies Mstar,CBOS./'^star,BCOi = 0.77 
and NEL-galaxies Mstar,CB0%iy^star,BC0Z = 0.70. This 
demonstrates that similar to the M05 fits, the differ- 
ences in Mgtar are not a strong function of SED type 
and are a systematic offset between the models. Of in- 
terest, this difference in Mgtar is less dramatic than the 
factor of 0.32 found by Bruzual (2007) when fitting the 
galaxies used in Maraston et al. (2006); however, they 
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Fig. 6. — SEDs fit with BC03 models (black), M05 models (green) and CB08 models (blue), plotted as vs. Log(A). Red points are 
the broadband photometric data and cyan points are the binned NIR spectroscopy. 



compared Mgtar using fits without including dust extinc- 
tion. Clearly, once dust is included in the fitting the 
differences in Mgtar are not as significant. 

The median r is similar between the CB08 and BC03 
fits; however, like the M05 fits, its dependence is strongly 
bimodal. In the CB08 fits, the EL-galaxies have r's 
which are 0.66 that of the BC03 r's. This behavior is 
similar to the M05 EL-galaxies, which also have shorter 
SF timescales. The NEL-galaxies actually have larger r's 



in the CB08 models compared to the BC03 models by a 
factor of 1.25. 

Unlike the M05 models, the (i)sFR of galaxies in the 
CB08 models are actually older than the BC03 models, 
with a median value 1.24 times larger. This dependence 
is again strongly bimodal, with the EL-galaxies having a 
(i)sFR 0.96 times the BC03 models, but the NEL-galaxies 
are older, with (t)sFR,CB08 = l-46(i)sFR,BC03- Interest- 
ingly, this gives the same relative behavior as the M05 
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Fig. 7. — As Figure 6. 

models, i.e., the EL-galaxies get younger relative to the 
NEL-galaxies; however, with the CB08 models it is the 
NEL-galaxies that get older while the EL-galaxies stay 
the same age. In the M05 models, the EL-galaxies get 
younger, but the NEL-galaxies stay the same age. 

The median extinction is lower in the CB08 models 
compared to the BC03 models, being AA„ — -0.4 for all 
galaxies. The effect is stronger for NEL-galaxies, which 
have AA„ = -0.6, compared to the EL-galaxies which 



have AA„ = -0.2. 

Overall, the net effect of using the CB08 models com- 
pared to the BC03 models is that the mean Mgtar of 
galaxies is lower by a factor of 0.74, and perhaps more 
importantly, the galaxy population becomes more bi- 
modal in several parameters. This is similar to the re- 
sult from the M05 models; however, the bimodalism is 
driven by movement in different types of galaxies. In 
the CB08 models compared to the BC03 models EL- 
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Fig. 8. — Left Panel: Comparison of Xr values for fits with and without the IRAC data using the BC03 models. Middle Panel: Same as 
left panel but for fits with the M05 models. Right Panel: Same as left panel but for fits with the CB08 models. Open red circles represent 
galaxies whose SEDs follow a power law in the IRAC bands of the form f,y oc ^ where ol < -0.5. These galaxies are commonly defined 
as "power-law galaxies" and likely have the rest-frame NIR, dominated by emission from dust heated by an AGN. The M05 models show 
the least increase in Xr when the rest-frame NIR data is included in the fitting; however, the BC03 and CB08 models have only a few 
non-power-law galaxies that have notably worse fits when the rest-frame NIR data is used in the fitting. 

6. SYSTEMATIC EFFECTS ON STELLAR POPULATION 
PARAMETERS FROM METALLICITY AND DUST LAW 
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Fig. 9.— Residuals from the SED fits using the BC03, M05, and 
CB08 models corrected to the rest-frame. The residuals are the 
fractional residual where A/a//a = fobs - f model/ fobs- The black 
line represents a running average of the nearest 30 data points. 
The running average suggests that the M05 and CB08 models may 
describe the SEDs of the galaxies better in the range 1.0 x 10^ 
A < A < 1.8 X 10^ A; however, the improvement is ~ 5-10%, well 
within the photometric errors. 



galaxies have lower r's, but similar (t)sFR 's, giving them 
a sHghtly burstier SFH. This is different than the M05 
models which are both significantly younger and have 
even smaller r's. In the CB08 models the NEL-galaxies 
have both longer r's, older ages, and surprisingly, less 
dust. In the M05 models, the NEL-galaxies had similar 
r's, and (i)sFR as the BC03 models. 

If we divide the galaxy population by SFR in both the 
M05 and CB08 model it is more bimodal than in the 
BC03 models. For the M05 models, the change is dom- 
inated by the EL-galaxies, which are notably younger 
and bustier with higher SFRs than their BC03 counter- 
parts. For the CB08 models, the bimodality is driven by 
the NEL-galaxies which become older, with longer r and 
smaller SFRs than their BC03 counterparts. 



As we have shown in the previous section, SED fits 
using the current generation of SPS codes provide stel- 
lar populations parameters that suffer from significant 
systematic effects. In this section we quantify how large 
these systematic effects are relative to other known sys- 
tematics such as the choice of dust law and metallicity. 

We allow three dust laws to redden the SEDS: the 
Milky Way (MW) and Large Magellanic Cloud (LMC) 
dust laws determined by Fitzpatrick et al. (1986), and 
the Small Magellanic Cloud (SMC) dust law determined 
by Prevot et al. (1984). We also incorporate two new 
metallicities into the fits, a subsolar (Z — 0.2 Zq), and su- 
persolar (Z = 2.5 Zq) population. For the BC03 models 
these are the m42, and m72 stellar populations, respec- 
tively. 

In an attempt to isolate variables we use the BC03 
models as our control SPS code. When comparing dust 
laws we assume solar metallicity as our control metal- 
licity, and when comparing metallicities we assume a 
Calzetti et al. (2000) dust law as the control dust law. It 
is possible that using a different control model will affect 
some of these comparisons; however, the large number of 
potential permutations allowed in choosing the control 
model (there are 36 possible combinations of SPS code, 
dust law, and metallicity) require us to chose only one 
in order to make the interpretation of systematic effects 
tractable. Moreover, the most commonly used models 
for determining stellar population parameters thus far 
have been the BC03 models with solar metallicity and 
the Calzetti dust law, so systematic changes in parame- 
ters relative to this model are probably the most inter- 
esting for understanding the implications of assumptions 
made in previous work. 

We first examine the effect of using the various 
dust laws on the parameter estimates. The ratios 
of parameters determined from fits using the SMC, 
MW, and LMC dust laws compared to those from the 
BC03_m62-|-Calzetti control model are listed in Table 5 
and are plotted in Figure 11. The open red triangles and 
open blue squares represent the medians of the NEL- 
galaxies and EL-galaxies, respectively, and the black 
points denote the medians of the whole sample. The 
error bars on the black points denote the range spanned 
by 50% of the sample. Figure 11 shows that when param- 
eters determined using the LMC and MW dust laws are 
compared to those determined using the Calzetti law, 
the only parameter that is significantly different is the 
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Fig. 10. — Ratio of best-fit CB08 (open triangles) and M05 (filled circles) SED parameters to the best fit BC03 SED parameters as a 
function of the BC03 pajrameters. The median M05 fit is 0.63 of Mst^r, and 0.65 of (t)sFR compared to the BC03 fits, while the z, SFR, 
and are similar. The median CB08 fit is 0.74 of Mgtari and 1-24 of (i)sFR! ''■'id has 0.4 mag less Av compared to the BC03 fits, while 
the 2, SFR are similar 
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Fig. 11. — Difference in the median value of parameters fit using the various SPS codes, dust laws, and metallicities (see text) compared 
to the "control" model. Open blue squares represent galaxies with detectable emission lines and open red triangles represent galaxies 
without detectable emission lines. The filled black circles are the median for all galaxies. The lines denote the range that contains 50% of 
the galaxies. Most of the largest systematic offsets are from the SPS codes, and their most significant effect is on Mgtai., and (t)gpi^, the 
two parameters best determined by the SED fitting. 
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A„. Parameters such as Mgtar, t, (t)sFR , and SFR only 
change systematically by of order ~ 10%, and therefore 
appear to be robust regardless of whether the MW, LMC, 
or Calzctti dust laws arc used. Furthermore, the changes 
do not appear to depend on the EL- and NEL-galaxy 
classification. As pointed out by Forster Schreiber et 
al. (2004), the lack of significant changes in parameters 
between these three dust laws is not surprising because 
they arc fairly similar in shape. The only notable differ- 
ence is that the Calzctti law lacks the 2175A absorption 
feature and that the ratio of total to selective absorption 
is slightly higher, = 4.05, compared to R„ = 3.1 for 
the MW and LMC laws. 

Using the SMC dust law docs cause some signifi- 
cant changes to several parameters. The median r and 
(i)sFR are larger by a factors of 2.5, and 2.8, respec- 
tively, and the median Ay is lower by 0.4 mag. The 
median SFR is similar; however, this results from a coin- 
cidental cancellation caused by the NEL-galaxics having 
SFRs a factor of ^ 2 higher, and EL-galaxies having 
SFRs a factor of - 2 lower. The SMC dust law has by 
far the strongest attenuation of all four dust laws at A < 
is 2000A rest-frame and therefore it requires stellar pop- 
ulations with more rest-frame UV flux to describe the 
SEDs. This leads to a preference for longer timescalcs 
for SF and older ages so that there is a contribution 
from both a young and old population simultaneously. 
Because of the shape of the SMC dust law, the fits prefer 
intrinsically redder stellar populations (i.e., those with 
older ages and longer r) and less dust, rather than short 
bursts with more dust. The preference for older ages 
and less dust would probably diminish if we used the 
SMC law with a subsolar metallicity population, which 
is intrinsically bluer and would allow for more dust and 
shorter timescales for SF. 

Next we examine the effect on the parameters caused 
by varying the metallicity while holding the SPS model 
and dust law fixed. The median values of the parame- 
ters determined from those fits compared to the control 
models fits are listed in Table 6 and are also plotted in 
Figure 11. The choice of metallicity affects more param- 
eters than does the choice of dust law; however, it is 
still a modest effect. The r of the stellar populations is 
lower for both metallicities; however, the dependence for 
EL- and NEL-galaxics is reversed. The A^, depends on 
metallicity in the predictable way, with subsolar metal- 
licity requiring more dust, and supersolar requiring less. 

Perhaps what is most interesting about the variations 
in parameters caused by choice of dust law and metallic- 
ity is that other than the SMC dust law, their net sys- 
tematic effect on Mstar, and (t)sFR , the two parameters 
that are the best-determined using SED fitting (§ 4.5) is 
fairly small, ~ 10-20%. They do have larger effects on 
T and SFR; however, these are the parameters that are 
most poorly determined using the SED fitting. 

One disconcerting aspect of Figure 1 1 is that it shows 
clearly that it is the SPS codes that have the largest 
systematic effects on most parameters, and most impor- 
tantly, that the largest systematic changes occur on the 
parameters that are best determined by the SED fitting: 
Mstar, (^)sFR j and Av. So, despite the fact that us- 
ing the best data currently available we can constrain 
the Mstar, (i)sFR , and A„ of individual galaxies to ~ 
0.1 dex, 0.3 dex, and ± 0.3 mag, respectively (§ 4.5), it 



appears that systematic errors in these parameters from 
the various SPS codes dominate over random errors for 
all three parameters as they cause systematic uncertain- 
ties of roughly ± 0.3 dex in Mstar and (f)sFR , and ± 
0.4 mag in A^. This suggests than unless these system- 
atic discrepancies between the SPS models are resolved, 
more multiwavelength data on the stellar component of 
the galaxy SED will not provide better constraints on the 
evolution of the stellar mass density in the universe, or 
the peak epoch of the formation of stars in the universe. 
Longer wavelength data that probes emission from dust 
(e.g., MIPS, ALMA) will almost certainly improve the 
SFR estimates; however, parameters determined from 
the stellar SED are now limited by systematic errors. 
Adding to the challenge is the analysis in § 5.2 that shows 
that even with the best constrained SEDs of young stel- 
lar populations currently available we will not be able to 
"fit out" these systematic problems, as all three models 
describe the SEDs equally well. 

7. IMPLICATIONS; CONFIRMATION OF THE EXISTENCE 
OF A SUBSTANTIAL POPULATION OF QUIESCENT 
GALAXIES AT Z ~ 2.3 

Using this sample of galaxies Kriek et al. (2006) have 
argued that ~ 50% of K-selccted galaxies at z ~ 2.3 
have strongly suppressed star formation based on their 
lack of observed emission lines and the fact that the rest- 
frame UV through optical SED was best described by 
an old quiescent population. Although they had the 
high-resolution NIR spectroscopy with which to model 
the stellar populations, Kriek et al. lacked the rest-frame 
NIR data from IRAC, without which it can be challeng- 
ing to distinguish between young-and-dusty and old-and- 
quiescent populations (e.g., Labbc et al. 2005; Papovich 
et al. 2006; Wuyts et al. 2007, § 4). Although unhkely, 
they could not completely rule out the possibility that 
some fraction of the systems without emission lines could 
be highly extincted star forming galaxies rather than old 
and quiescent galaxies, and that the lack of emission lines 
and red optical colors are the result of significant dust ob- 
scuration. Here we examine the stellar populations from 
the SED fits of the EL and NEL galaxies to verify that 
the SEDs of the galaxies without emission lines are still 
consistent with them being systems with strongly sup- 
pressed star formation. 

In the top panels Figure 12 we plot the specific star 
formation rate (SSFR), defined as SFR/Mgtar, as a func- 
tion of Mstar for the EL galaxies (blue) and NEL galaxies 
(red). Candidate AGNs selected using the emission-line 
diagnostics in Kriek et al. (2007) as well as PLCs (§ 5.2) 
are indicated with an open diamond. As we showed in § 
6, the SPS codes have the largest systematic effects on 
Mstar and SFR, and therefore we plot three versions of 
the diagram in Figure 12, one each for values determined 
using the three SPS codes. 

It is clear from Figure 12 that the majority of galaxies 
without emission lines are found at the lowest SSFRs, 
whereas those with emission lines are found at the high- 
est SSFRs. There are some notable mo del- dependent dif- 
ferences, with the M05 and CB08 models implying that 
the population of galaxies at z ^ 2.3 has a more bimodal 
distribution of SSFRs, whereas the BC03 models suggest 
more of a continuum in SSFRs. This is apparent in the 
bottom panels of Figure 12, where we plot histograms of 
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the SSFRs. 

If we adopt the Kriek et al. (2008b) definition of qui- 
escent systems as those that have SSFRs < 0.05 Gyr~^ 
(i.e., those that will increase their Mgtar by < 5% in the 
next Gyr if the SFR remains at the same value) , then we 
find that 13, 17, and 22 of the 34 galaxies (38%±12%, 
50%±15%, 65%±18%) would be classified as quiescent 
based on the SED parameters from the BC03, M05, and 
CB08 models, respectively. These numbers arc in good 
agreement with the 45^12% determined by Kriek et al. 
(2006) based on the fraction of galaxies without emis- 
sion lines and SED fitting. Despite the good agreement, 
there are a few subtle issues in the emission line classi- 
fication. There are 16 galaxies without emission lines, 
and of these, 10, 12, and 15 of these would be considered 
quiescent based on their SED parameters from fits to 
the BC03, M05, and CB08 models, respectively. There- 
fore, the majority of systems without emission lines are 
quiescent; however, there is a non-negligible population 
(~ 6% - 33%) of galaxies without emission lines that do 
have some ongoing obscured star formation ranging from 
^ 10-850 Mq yr~^. Coincidentally, these few galaxies are 
canceled out in the overall quiescent fraction by a small 
fraction of emission line galaxies that are likely to be qui- 
escent based on their SED, but probably have an AGN 
component that is responsible for the emission lines. 

Taken together, these results suggest that the classifi- 
cation of 2: ~ 2 galaxies into star forming and quiescent 
categories based on the presence/lack of emission lines 
and without any information on their SED is probably 
correct for the majority of systems (~ 70-90%); however, 
there is some cross contamination between the two cate- 
gories from dusty star forming galaxies without emission 
lines, and quiescent galaxies with an emission line AGN. 
Within our sample these two populations appear to be 
small, and equally abundant and therefore cancel each 
other when estimating the quiescent fraction. As a re- 
sult, the fraction of quiescent galaxies is roughly the same 
using the SED fits or the emission line classification. Our 
sample is fairly small, so it is not clear if the same can- 
cellation would occur for larger samples. 

Although the contamination in the emission line clas- 
sification suggests we should appeal to the SED fits for a 
more robust classification, that method is not completely 
"clean" cither: the quiescent fraction ranges by roughly 
a factor of 2 depending on which SPS code is used to 
model the SEDs. Considering this systematic effect and 
using the la confidence intervals from the models, we 
can conclude that the fraction of K-selectcd galaxies at 
2: ~ 2.3 that are classified as quiescent based on their 
SEDs ranges from as low as 26% to as high as 83% at 
68% confidence, in good agreement with the value de- 
termined from the emission line fraction by Kriek et al. 
(2006). Clearly this uncertainty is large; however, even 
within our modest sample of 34 galaxies, the statistical 
uncertainties are already smaller than the systematic un- 
certainty from the SPS codes. 

The histograms in Figure 12 suggest a more bimodal 
distribution of SSFRs when using the M05 and CB08 
models compared to the BC03 models. Kriek et al. 2008b 
have shown that there is already a developed red se- 
quence in the (U-B)rest vs. Mgtar plane for this sample of 
galaxies and that most of the galaxies on the red sequence 



have low SSFRs as inferred from SED fits using the BC03 
models. It is enticing to prefer SPS models that show the 
same bimodality in SSFRs as in rest-frame color; how- 
ever, given that the current sample is small, and the ran- 
dom uncertainties in the SSFRs are large (factor of 2-3) 
it is better to defer such a comparison to a larger sample 
of galaxies (e.g., van Dokkum et al. 2009). It will be in- 
teresting to see how the SSFRs of galaxies in the larger 
sample relate to their rest-frame colors, which will be a 
much more accurately determined quantity using those 
data. 

8. SUMMARY 

In this paper we performed SED fitting on a sample 
of 34 K-sclccted galaxies at 2 ~ 2.3. With NIR spec- 
troscopy as well as deep photometry in thirteen broad- 
band filters, these galaxies have arguably the best con- 
strained SEDs of 2; > 2 galaxies currently available. With 
this sample we studied the importance of the rest-frame 
NIR in determining stellar population parameters, the 
systematic differences between the current generation of 
SPS codes, as well as quantify the systematic eflfects 
in determining stellar population parameters of young 
galaxies from the choice of SPS code, dust law, and 
metallicity. Our main results are summarized below: 

1. Using the BC03 models with solar metallicity, the 
Calzetti dust law, and a Salpeter IMF as a control model, 
fits with and without the IRAC data show that includ- 
ing the rest-frame NIR in the SED fitting of young stellar 
populations provides improved constraints on nearly all 
stellar population parameters. In particular, r, A^,, and 
SFR are improved by factors of ~ 1.5 - 2.0, and Mgtar is 
improved by a factor of 1.3. 

2. The IRAC data is particularly helpful when only 
broadband data (no spectroscopy) is available at other 
wavelengths. In particular, constraints on Mgtar, SFR, 
and A„ improve by factors of 4, 2.5, and 0.5 mag, re- 
spectively. 

3. IRAC data only modestly improves the photometric 
redshifts of 1.6 < z < 2.9 galaxies. The improvement 
only occurs for galaxies with low S/N JHK photometry, 
and only when all four IRAC bands are included in fit- 
ting. Our comparisons show that deep JHK data is far 
more valuable than IRAC data when determining the 
2;phot of galaxies at 1.6 < z < 2.9. 

4. Neglecting systematic effects from choice of SPS code, 
dust law, metallicity, and IMF, and using our entire pho- 
tometric dataset (U— >8/im-|-NIRspec), we find that the 
best possible constraints on Mgtar, {t)sFR, and A„ of in- 
dividual galaxies in our sample from SED fitting are ± 
0.12 dex, ± 0.26 dex, and ±0.3 mag, respectively. Even 
with these well-constrained SEDs, the r and SFR of in- 
dividual galaxies are still uncertain to factors of 2-3. Us- 
ing only the broadband data (U^8/um) as well as Zphot, 
the constraints on Mgtar , (^)sfr, and A„ are ±0.11 dex, 
±0.36 dex, and ±0.3 mag, nearly as good as when the 
NIR spectroscopy is included in the fitting. 

5. Comparison of the quality of fits provided by the 
BC03, M05, and CB08 models using several methods 
shows that all three provide equally good descriptions of 
the data. There is tenuous evidence that the BC03 and 
CB08 models describe the rest-frame UV through optical 
SED slightly better than the M05 models; however, the 
latter models may provide a slightly better description of 
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Fig. 12. — Top Panels: Specific star formation rate (SSFR = SFR/Mgtar) ^ ^ function of Mgtar determined from fits to the BC03, M05, 
and CB08 models. Galaxies with detected emission lines are plotted in blue, those without detectable emission lines are plotted in red. 
Candidate AGN (see text) are indicated using open diamonds. The three SPS models produce significantly different distributions in this 
parameter space; however, galaxies without emission lines typically have much lower SSFRs than those with emission lines. Bottom Panels: 
Histograms of the SSFR determined using models from the three SPS codes. The dashed line represents our definition of a quiescent galaxy 
(SSFR < 0.05 Gyr-i). Based on this definition, 38%±12%, 50%±15%, 65%±18% of our sample would be consider quiescent using the 
BC03, M05, and CB08 models, respectively 



the rest-frame NIR SED. Taken together, there is httle 
evidence that one set of models provides an overall better 
description of the SEDs of young massive galaxies based 
on the quality of fit. Furthermore, given that the differ- 
ences in the predicted fluxes of the best-fit models at all 
wavelengths is <5%, it is unlikely that the current gener- 
ation of data on the SEDs of young, massive galaxies can 
inform the models enough to aid in their improvement. 

6. Comparison of the systematic changes in stellar popu- 
lation parameters determined using the three SPS codes 
shows that when we fit with the M05 models the Mgtar 
and (t)sFR of galaxies go down by factors of 0.63 and 
0.65, respectively relative to the BC03 fits. Galaxies fit 
with the CB08 models have 0.75 the M^tar of the BC03 
fits, but are on average older by a factor of 1.24. No- 
tably, the changes in most parameters appear to depend 
on whether galaxies have detectable emission lines. Be- 
cause they affect the SFRs rates of SF and NEL-galaxies 
differentially, both the M05 and CB08 models provide an 
overall stellar population that has "burstier" episodes of 
star formation than the BC03 models, as well as suggest 
a galaxy population that is more bimodal in terms of the 
overaU SSFRs. 

7. Comparison of systematic changes in stellar popula- 
tion parameters caused by choice of dust law, metallicity, 
and SPS code show that the latter are the most signifi- 
cant of the three. The systematic effects on Mgtar, and 



(^)sFR, the best determined parameters from the SED 
fitting, from different dust laws and metallicities are gen- 
erally small, of the order of 10-20% on Mgtar, and (i)sFR- 
Conversely, the systematic effects from the SPS code on 
these parameters are most significant. These systematic 
effects are larger than the random errors from the SED 
fitting and suggest that further progress in determining 
the evolution of the stellar mass density and the peak 
epoch of SF in the universe will require significant im- 
provements in the models. 

8. Comparison of the SSFRs determined using the SED 
fitting for galaxies with and without detectable emission 
lines shows that most galaxies without emission lines are 
quiescent (66% to 96% depending on SPS model); how- 
ever, the remainder are likely to star forming galaxies 
that have obscured emission lines. In our sample, the 
obscured emission line population is roughly equal to the 
population of galaxies that are quiescent based on their 
SED, but have detected emission lines, probably from 
an AGN. Therefore, despite some small cross contamina- 
tion of quiescent and active galaxies in the emission line 
and non-emission line categories, respectively, the over- 
all fraction of K-bright, z ~ 2.3 galaxies that would be 
considered quiescent based on their SED is ~ 50% (26% 
- 83% depending on SPS model), nearly identical to the 
fraction without emission lines. 
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APPENDIX 

DETERMINATION OF SYSTEMATIC AND RANDOM ERRORS IN STELLAR POPULATION PARAMETERS 

In Figure 13 we plot histograms that show the distributions of parameters determined with the U^K, 
U^z'-|-NIRspec, and U— +8/im data relative to those determined with the U^8/im-|-NIRspec data. The mean and 
median of each distribution is listed in the upper right of each panel. These demonstrate the systematic differences in 
the parameters with various datasets relative to the U— s-8/im+NIRspec data, and are the parameters listed in Table 1. 

The standard deviation (cr) of the distributions are listed in the upper right of each panel in Figure 13. These 
serve as an indication of the scatter between the best fit parameters, even though some of the distributions are non- 
Gaussian. Standard deviations are sensitive to large outliers and therefore we also compute the normalized median 
absolute deviation (ctnmad) to compare if the a is representative of the true scatter in parameters or is dominated by 
a few large outliers. For a Gaussian error distribution, cnmad = rms. 

Because the histograms are ratios of parameters computed with different datasets, the scatter is caused by the 
combination of errors in both datasets. To determine whether the given data improves the quality of the fits we must 
assess how much of the scatter is caused by uncertainties in U^K, U^z'-f NIRspec, or U^8/xm fits and how much is 
intrinsic to the SED fitting with the comparison sample (U— s-8/xm-l-NIRspec). The intrinsic contribution is estimated 
by simulating a new distribution of the U^8/^m+NIRspec fit parameters based on the error estimates listed in Table 
1. We generate this distribution by drawing 200 new values of each parameter, for each galaxy, using the uncertainties 
and compare the new parameters to the original ones. The resulting "expected" distributions are plotted in blue in 
Figure 13 with the measured mean, median, cr and unmad also listed in blue. These values are the scatters used on the 
X-axis of Figure 13. Where the grey and blue histograms differ in mean, median, or sigma, the IRAC data provides 
additional constraints. The main results from the histograms in Figure 13 are summarized in § 4.3 and § 4.4. 
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Fig. 13. — Left Panels: Comparison of stellar population parameters determined using only the broadband U— »K photometry to those 
determined using all of the data (U— >8/im-f NIRspec). Middle Panels: As the left panels, but comparing parameters from broadband plus 
NIR spectroscopy (U— >z'+NIRspec) to those determined with U— >8/.tm-f NIRspec. Right Panels: As left panels but comparing parameters 
from broadband optical and IRAC photometry (U— >8/xm) to those determined with U— >8/.tm+NIRspec. The blue curves in each panel 
are the spread in the distribution that would be expected based on the Monte Carlo errors determined using the U^8/im+NIRspec 
data. If the parameters determined from a subset of the data have the same accuracy and precision as those determined from the entire 
U^8/.tm+NIRspec data set, the blue and grey histograms should be similar. The differences in the accuracy of the parameters using the 
various datasets can be determined by comparing the a or the tj^jj^j^q accross the different datasets. The addition of rest-frame NIR 
data from IRAC provides significant improvement in most parameters, particularly compared to those determined from U— >K photometry 
alone. The IRAC data provides better constraints than the NIR spectroscopy on the parameters Mst^n Aj,, and SFR; whereas the NIR 
spectroscopy provides better constraints on r and (t)sFR- 
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TABLE 1 

SED Modeling Parameters 



Name 
(1) 



Model 



(2) 



(3) 



age 
Gyr 
(4) 



Gyr 
(5) 



mag 
(6) 



M, 
10" M, 
(7) 







SFR 

Mq yr-i 
(8) 



< tSFR > 
Gyr 
(9) 



<Xr> 



1030-32 



BC03 (no irac) 
BC03 

MAR05 (no irac) 
MAR05 
CB08 (no irac) 
CB08 



2.23 
2.23 
2.23 
2.23 
2.36 
2.36 



+0.13 
0.02 
+0.13 
0.03 
+0.13 
0.03 
+0.02 
0.14 
+0.00 
0.14 
+0.00 
0.14 



1.00 
0.90 
1.00 
0.90 
0.90 
0.90 



+U.OV 
0.12 
+0.07 
0.09 
+0.20 
0.06 
+0.08 
0.09 
+0.12 
0.04 
+0.12 
0.04 



0.00 
0.20 
+0.60 
0.20 
+0.25 
0.20 
+2.10 
0.10 
+0.10 
0.20 
+0.10 
0.00 



0.01 
0.03 
0.05 
0.01 
0.01 
0.01 



0.10 

+0.40 
0.10 
1.75 
0.80 
+0.30 
0.40 
+2.60 
0.50 
1.45 
0.30 



+0.09 
0.00 
+0.12 
0.02 
+0.10 
0.04 
+0.39 
0.00 
+0.09 
0.00 
+0.07 
0.00 



0.1 
0.4 
0.0 
0.0 
0.1 
0.0 



+0.27 
0.01 
+0.08 
0.03 
+0.28 
0.11 
+0.11 
0.03 
+0.38 
0.11 
+0.20 
0.09 



+0.4 
0.1 
+0.3 
0.4 
+0.1 
0.0 
+0.0 
0.0 
+0.4 
0.1 
+0.0 
0.0 



1.34 
1.63 
1.16 
1.05 
1.46 
1.30 



+0.3 
-1.3 
+0.2 
0.4 
+ 1.7 
0.4 
1.1 
0.3 
+0.7 
1.0 
+0.3 
0.9 



+0.30 
0.24 
+0.25 
0.25 
+0.12 
0.20 
1.44 
0.13 
+0.41 
0.24 
+0.11 
0.05 



+0.B2 
0.54 
+0.54 
0.49 
+1.22 
0.44 
+0.81 
0.18 
+0.93 
0.70 
+0.45 
0.82 



hO.O 
-0.0 
-1+0.5 
-0.0 
,+0.2 
-9.6 

n 0+''-'' 
"J-"-o.o 

-,+0.0 
-0.0 
-1+0.0 
-0.0 



o.o: 
o.oj 

9.6^ 



0.01 

o.oj 



0.99 
0.87 
0.95 
0.89 
0.89 
0.89 



n,+i3.b 
'^-O.O 

"^•"^-0.0 

-1.2 
+6.6 
-0.9 



0.00 
0.21 
+0.60 
0.20 
+0.15 
0.18 
1.71 
0.13 
+0.10 
0.20 
+0.10 
0.04 



n 07+1.3::! 

-0.09 
36+016 

0.88^ ' 
0.66 
0.78 
0.70 



1.56 
1.46 
1.92 
1.89 
1.55 
1.48 



1030-101 



BC03 (no irac) 
BC03 

MAR05 (no irac) 
MAR05 
CB08 (no irac) 
CB08 



1.79 
1.79 
1.72' 
1.79 
1.72 
1.72 



0.40 
0.40 
1.00 
0.70 
0.90 
0.80 



0.02 
0.04 
0.12 
0.04 
0.12 
0.10 



1. 
1.9 
0.6 
0.7; 
1.1 
1.2 



2.76 
2.91 
1.45 
1.49 
2.10 
2.14 



O.O"* 



1-2Z 
0.9 



T-Bir 

0.85 
+ 1.35 
0.70 
+0.30 
0.55 
+0.40 
0.05 
+0.70 
0.60 
0.00 
0.04 



J^'J.U_g gg 

10 o+i° o 

-1+0.00 
-19.9 
,+19.8 
-0.09 
5 00+15.0 

1+0.01 

-0.03 



-+U.4 
-0.3 
-+0.3 
-0.3 
^+0.6 
-0.2 
a+0.1 
-0.4 

1.5+0-3 



0.7010;°^ 

1+0.23 

-0.08 



ii8-2:iS5;2 
-100.6 
-53.1 

140.61396^2 

onq O + 126.0 

OUO.O_j^gg 3 

118 4+104.2 

131.8tj^g'^g 



1.47 
-0.70 
+0.33 
0.38 
+2.22 
0.43 
1.25 
0.22 



45+o'«i 

45+0-''i 

30+0- 16 

05+0-20 

41+0-41 
"■^1-0.28 

06+0-00 



0.74 
0.82 
0.84 
0.89 
0.80 
0.81 



1030-609 



BC03 (no irac) 
BC03 

MAR05 (no irac) 
MAR05 
CB08 (no irac) 
CB08 



1.800 
1.800 
1.800 
1.800 
1.800 
1.800 



0.90 
0.90 
0.60 
0.10 
0.80 
0.10 



20.0: 
0.12"* 



0.04'' 



1.5"' 
I.5-* 
1.5^ 
1.9^ 



0.41"' 



118.2"* 



1, 



TTTT 
0.20 
0.30 
0.30 
+0.10 
0.54 
+0.10 
0.51 
+0.20 
0.20 
0.20 
-0.10 



12+0-03 
^•1^-0.04 
-,+0.05 
-0.09 



-0.3 
a+0.0 
-0.2 



0.82iO-J4 



0.50 



+0.3 
0.2 
+0.5 
0.4 
1.8 
0.0 
+1.5 
0.0 
+0.6 
0.2 
+0.0 
0.2 



+0.01 
0.09 



+0.17 
0.12 
+0.22 
0.25 
+0.30 
1.52 
+0.30 
0.00 
+0.34 
0.14 
+0.00 
0.13 



2.4 
+ 14.3 
12.9 
+28.4 
0.1 
+2.1 
0.1 
+9.1 
2.7 
+0.0 
2.6 



+u.r/ 

0.08 
+0.25 
0.23 
+0.06 
0.47 
+0.06 
0.44 



1.09 
1.12 
1.12 
1.19 
1.02 
1.23 



1030-807 



BC03 (no irac) 
BC03 

MAR05 (no irac) 
MAR05 
CB08 (no irac) 
CB08 



2.367 
2.367 
2.367 
2.367 
2.367 
2.367 



0.70 
0.50 
0.60 
0.60 
0.70 
0.70 



0.10"* 



0.08 
0.08 
0.12 
0.12 



0.00 
+0.00 
0.75 

+0.60 

0.10 
+0.10 
0.20 

+0.50 
0.00 
+0.50 
0.30 



+0.04 
0.07 
+0.04 
0.07 
+0.03 
0.04 
+0.03 
0.04 



0.2 
0.8 
0.0 
0.0 
0.2 
0.2 



1.11 
1.47 
0.82 
0.82 
1.11 
1.10 



3.6 
12.9 
0.7 
0.7 
3.6 
3.6 



0.58 
0.40: 
0.52 
0.52 



0.581q;i7 



0.58 



+0.80 
0.00 
+10.0 
7.50 
+9.75 
0.15 
17.5 
1.85 
+0.80 
0.00 
1.60 
0.50 



+0.2 
0.1 
+0.1 
0.0 
+0.1 
0.2 
+0.0 
0.1 
+0.2 
0.1 
+0.0 
0.1 



1 43+U.VO 
l-^'^-O.OO 
o 20+0.40 

-1+0.57 



-1+69.7 
-28.4 
-+38. 

-3.7 
7+61.4 
-90.2 

219-lli8'6 

+69.7 
28.5 
+4.3 
33.3 



+0.17 
0.06 



34+0.24 
'J-'^^-O.OO 
-,+0.05 
-0.33 

0^25+0-25 



1.87 
1.88 
1.97 
1.78 
1.87 
1.74 



1030-1531 



BC03 (no irac) 
BC03 

MAR05 (no irac) 
MAR05 
CB08 (no irac) 
CB08 



2.613 
2.613 
2.613 
2.613 
2.613 
2.613 



0.50 
2.50: 
0.40 
0.90: 
0.50 
1.00 



0.20 
10.0 
0.25 
2.50 
0.20 
0.90 



0.4 
0.6 
0.6 
0.7; 
0.4 
0.6 



1.30 
1.91 
1.43 
2.14 



+0.14 
0.15 
+0.04 
0.25 
+0.14 
0.13 
+0.14 
0.25 
+0.14 
0.14 
+0.27 
0.00 



0.30 
2.60 
0.40 
+0.80 
0.65 
+3.20 
0.22 
1.65 
0.40 
0.00 
2.50 



+0.28 
0.11 
+0.60 
0.12 
+0.18 
0.11 
+0.87 
0.02 
+0.28 
0.11 
+0.05 
0.45 



+0.5 
-1.1 
+0.4 
0.9 
+1.9 
0.5 
+0.8 
1.3 
+0.7 
1 

+0.5 
0.3 



0.21 
+0.17 
0.23 
+0.70 
0.00 
+0.27 
0.40 



81.03 
153.5j 
160.2j 



81.0 
149.7 



1.30^ 



0.47; 
0.34; 
0.59 



1 n2+0-E>V 
l-y^-0.55 

1.72+1-34 



+81.9 
32.4 
+8.7 
-10.6 
-123.6 
-3.1 
+67.3 

0.2 
+60.9 
-12.2 
+2.6 
0.9 



0.03 
+0.05 
0.09 
+0.24 
0.00 
+0.24 
0.16 



+0.98 
0.21 
+2.15 
-0.30 
+0.67 
0.51 
+2.47 
0.19 
+0.72 
0.23 
+0.04 
2.05 



1.38 
1.81 
1.58 
1.60 
1.39 
1.36 



1030-1611 



BC03 (no irac) 
BC03 

MAR05 (no irac) 
MAR05 
CB08 (no irac) 
CB08 



1.93 
1.93 
1.93 
1.93 
1.93 
1.69 



0.50 
0.90 
0.70 
0.30: 
0.60 
3.50 



0.12 
0.20 
0.12 
0.03 
0.12 
0.65 



1 

1.2 

0.7; 
1.4 
1.5 
0.3 



0.98 
1.19 
1.63 
1.62 



T73ir 
1 .55 
0.20 
0.40 
1.25 
1.10 
+0.10 
0.10 
+0.50 
1.65 
+0.20 
0.10 



+ 18.0 
1.75 
+0.15 
0.15 
+ 17.2 
2.63 
+0.04 
0.03 
18.0 
1.80 
+0.08 
0.02 



+0.3 
0.3 
+0.3 
0.2 
+0.3 
0.5 
+0.3 
0.7 
+0.3 
0.3 
+0.2 
0.3 



0.39 
+0.82 
0.18 
+ 1.67 
0.43 
+0.58 
0.41 
+0.32 
0.68 



32.4 
12.1 
3.1 
0.2 
12.2 
1.5 



0.38 
0.71 
0.58; 
0.27; 
0.48 



+0.77 
1.36 
+0.19 
0.35 
+0.91 
1.05 
+0.00 
0.36 
+0.78 
1.40 
+0.30 
0.17 



+ 133.7 
-100.6 
+88.5 
34.7 
+94.5 
-51.2 
+ 14.8 
4.3 



-, qo+U.44 

1.00_Q gg 

0-581028 

0.8ll»j« 

99+0.09 
^^-^^-0.08 

1 qO+0.44 

n /I8+013 
0-48„o 09 



0.75 
0.81 
0.75 
0.73 
0.77 
0.96 



1030-1839 



BC03 (no irac) 
BC03 

MAR05 (no irac) 
MAR05 
CB08 (no irac) 
CB08 



2.312 
2.312 
2.312 
2.312 
2.312 
2.312 



2.25 
0.80 
1.50 
0.30 
2.25 
0.60 



2.00 
0.25 
2.75 
0.08 
2.00 
0.12 



1.1 
0.7: 
1.2 
1.0 
1.1 
0.3 



2.94 
1.50 
2.14 
1.15 
2.94 
0.92 



136.4 
43.0 
93 



0.30 
+0.20 
0.30 
+0.00 
0.20 
+0.10 
0.10 
0.50 
0.30 
+0.20 
0.10 



+0.05 
0.05 
+0.05 
0.08 
+0.00 
0.03 
+0.00 
0.03 
+0.05 
0.05 
+0.05 
-0.00 



7+0-4 
'^■'-0.7 
,+0.3 
-0.2 
,+0.3 
-0.1 



+0.50 
1.07 
+0.46 
0.41 
+0.31 
0.25 
+6.10 
0.27 
+0.50 
0.80 
+0.39 
3.0 



+ 14.2 
11.0 
+33.0 
19.1 
+4.9 
1.7 
+0.0 
-4.8 
12.1 
9.9 
11.2 
0.0 



0.25 



2.74 
2.73 
2.71 
2.74 
2.74 
3.01 



1030-2026 



BC03 (no irac) 
BC03 

MAR05 (no irac) 
MAR05 
CB08 (no irac) 
CB08 



2.511 
2.511 
2.511 
2.511 
2.511 
2.511 



1.00 
0.80 
0.90: 
0.80 
1.00 
0.80 



0.20 
0.20 
0.15 
0.15 
0.20 
0.15 



1.2: 
0.2"* 



0.5 
0.7; 

0.7; 



(-0.2U 



+U.Ob 
-0.14 
+0.05 
0.14 



+0.0 
0.3 
+0.3 
0.5 
+0.2 
0.1 



3.01 
3.70 
1.80 
2.10 
2.99 
2.55 



13.6 

44.8 
3.8 
8.7' 

13.5 

10.8 



0.80" 



61+0-15 

U.Di_o 22 



0.75 
0.65 
0.80 



+0.00 

0.17 

+0.09 
0.07 
+0.44 
0.25 



7+1-^ 
'-'-'-0.6 
7+0.6 

-0.1 



1 54+U.38 
1-^^-0.40 
-+0.23 

-0.24 



r+32.0 

^-6.5 
+8.3 
6.5 



65+0-15 
'J-PO-0.09 



n65+oiy 

'J-D3_0.49 
-+0.15 
-0.38 



1.46 
1.77 
1.58 
1.81 
1.44 
1.66 



1030-2329 



BC03 (no irac) 
BC03 



2.236 
2.236 



0.80^ 



0.15 



0.7"' 



1.55"' 



6.5 



0.65^ 



1.38 
1.29 



25 



TABLE 1 



Continued 



Name 
(1) 



Model 
(2) 



(3) 



age 
Gyr 
(4) 



Gyr 
(5) 



mag 
(6) 



M* 
10"Mq 
(7) 



SFR 
Mq yr- 
(8) 



< tSFR > 

Gyr 
(9) 



<Xr> 



MAR05 (no irac) 
MAR05 
CB08 (no irac) 
CB08 



2.236 

2.236 

2.236 

2.236 
r+u.ia 

^-0.04 
^+0.19 
-0.04 



0.70 
0.70 
0.70 
0.80 



-0.10 
-0.10 
-0.20 
-0.30 
-0.50 
0.20 
-0.10 



9n+»-«u 

"•'^'^-0.50 

-0.80 
-0.20 
-0.60 
-0.30 
-0.60 
-0.20 
-0.70 
-"-(). 30 

1 00+"-"" 
^■^'^-0.30 



U.1U_Q 02 
'J- '-'-0.09 
'J-J^^-0.11 

r, ir,-|-0.03 

'-'•^^-0.02 



0.1 

0.1 
0.6 
0.3 



-1-0.08 
0.11 
-1-0.12 
0.07 
-1-0.10 
0.09 
+0.10 
0.09 
+0.10 
0.09 
+0.00 
0.11 



+0.3 
0.1 
+0.2 
0.1 
+ 1.2 
0.3 
+0.2 
0.3 



0.86 
0.85 
1.32 
1.12 



0.6 
+0.4 
0.6 
+0.5 
0.0 
+0.2 
0.0 
+0.6 
0.6 
+0.3 
0.2 



+0.10 
0.08 
+0.05 
0.11 
+0.59 
0.19 
+0.16 
0.10 



1 0+"'2 

^^•'-'-0.2 
-1+0.2 
-1.0 
4^3+27.9 



l.Ot 



1.5 



+O.VU 
0.51 
+0.37 
0.37 
+0.40 
0.13 
+0.43 
0.15 
+0.80 
0.40 
+0.14 
0.33 



4.3 
+2.2 
0.3 



0.60 
0.60 
0.58; 
0.68 



0.11°;? 
0.lt°;« 



-0.08 
0.08 
0.14 
+0.22 
0.41 
+0.17 
0.09 



+U.b2 
0.30 



1.25 
1.13 
1.38 
1.33 



1030-2559 



BC03 (no irac) 

BC03 

MAR05 (no irac) 
MAR05 
CB08 (no irac) 
CB08 



2.25Z 
2.25j 
2.30"* 



-,+0.08 
-'-0.08 
2 of,+0.14 

1+0.14 
■'-0.06 

^••^"-0.06 



2.30"* 



0.70: 
0.90l 
O.QOl 
0.80^ 



0.12 
0.08 
0.10 
0.10 
0.10 
0.12 



0.6 
0.8 
0.0 
0.0 
0.6 
0.2 



1.84 
1.92 
1.13 
1.12 
1.83 
1.49 



0.78 



62+° ®'' 



-0.10 
-0.20 
+0.20 
0.00 
+0.30 
-0.10 
-0.10 
-0.00 
-0.10 
-"-0.20 

'^■o'^-o.io 



+0.02 
-0.07 
+0.06 
0.01 
+0.05 
0.01 
+0.03 
0.04 
+0.02 
-0.07 
+0.02 
0.06 



+0.4 
-0.3 
+0.1 
0.3 
+0.5 
0.6 
+0.1 
0.2 
+0.4 
-0.3 
+0.1 
0.1 



+0.47 
0.34 
+0.27 
-0.11 

1.70+°-*8 



'^■^-6.1 
7+29.1 
^•'-0.7 
0.0 

1.2 



-1+0.50 
-0.22 

0.80l°:50 
0.70l°-«° 
0.881°;°^ 



0.80^ 



"-0.15 
1O+0.14 

0^27+0.24 



1.39 
1.26 
1.57 
1.47 
1.39 
1.32 



1030-2728 



BC03 (no irac) 
BC03 

MAR05 (no irac) 

MAR05 
CB08 (no irac) 
CB08 



2.504 
2.504 
2.504 
2.504 
2.504 
2.504 



0.50 
0.20: 
0.30 

0.40 

o.5o: 



0.08 
0.02 
0.02 

0.05 
0.08 
0.08 



0.8 
1.6 
0.9 

0.6 
0.8 
0.8 



1.93 
2.48 



0.42: 



1.49 
1.93 
1.92 



+0.14 

0.13 
+0.11 

0.14 
+0.09 
0.07 
+0.12 
0.08 
+0.14 
0.13 
+0.12 
0.14 



,-1-0.85 
0.31 
,+0.60 
0.30 

+0.10 
0.15 
0.20 
-0.05 
0.90 
0.01 
+0.90 
0.01 



.,+0.04 
-0.03 
,+0.07 
'^-O.Ol 
1+0.24 
-0.00 
1+0.24 
-0.00 



+0.4 
1.3 

^•*-0.7 
+0.6 

0.4 
+0.2 
-0.5 
1 7+0.2 
-1.5 
+0.0 
1.5 



0.51 
+0.15 
0.19 
+0.22 
-0.34 
+0.08 
0.32 



+10.9 
0.0 

+3.8 
1.2 



6.1 
6.l^ 



+0.7 



0.35 
0.42 
0.42 



+0.3Y 
0.28 
+0.37 
0.26 
+0.12 
0.29 
+0.22 
-0.23 
+0.32 
0.07 
+0.25 
0.04 



^^■"-41.0 

41 n+30-4 

+ 77.8 
17.0 
+ 58.7 
64.7 
+39.6 

0.0 
+ 11.3 
0.0 



0.09 

+0.07 
0.00 
0.09 

-0.14 
0.08 
0.07 



,+O.Vi 
0.20 
,+0.47 
0.20 

+0.08 
0.10 
0.15 
-0.03 
+0.71 
0.00 
+0.71 
0.00 



1.93 
1.94 
2.04 

1.96 
1.93 
1.76 



1030-2927 



BC03 (no irac) 
BC03 

MAR05 (no irac) 
MAR05 
CB08 (no irac) 
CB08 



1.83: 
1.83 
1.84' 
1.83 
1.83 
1.83 



0.40: 

o.4o: 

0.20 
0.10 
0.10 
0.10 



0.04: 

0.02^ 
O.Olj 



1.4-f 



1.3 
1.8 



1.7"f 



1.05 
1.05 

0.83 
0.74 
0.79 
0.79 



17.3 
65.8 
0.4 
0.4 



0.29; 
0.29; 

0.16 
0.07 
0.09 
0.09 



3+U.2U 
-0.18 
3+0.15 
-0.21 
-1+0.18 
-0.03 
.+0.18 
■'-0.03 
9 nq+0.20 
^■^^-0.03 
3+0.20 
-0.20 



3+1.30 
-0.00 

-1+1.30 
-0.40 
.,+0.30 
-0.61 
-0.35 
-0.80 

1+1.76 
-0.40 



1+0.24 

-0.00 
,+0.18 
--0.11 
-1+0.05 
-0.09 
,+0.08 
-0.11 

U-1U_Q gg 

-+0.25 
-0.14 



+0.3 

1.8 
+0.5 
-0.8 
+ 1.7 
0.4 
+ 1.9 
0.1 
1.0 



*-l^-1.51 

'-'•"^-0.70 
5+1.85 

-0.55 

3.60ti;33 

,+1.94 
-0.66 



+ 194.9 
0.0 
+23.4 
11.9 
+ 19.5 

2.6 
+ 16.6 
-1.1 
"' 3 



+3.1 
3.1 



a+l-ii 
-0.02 

n + 0.87 

'-0.31 
-1+0.27 
-0.41 
1-0.09 
-0.69 

0.50l°;li 

-+1.60 
-0.25 



1.20 

1.14 

1.06 
1.08 
1.32 
1.23 



1256-0 



BC03 (no irac) 

BC03 

MAR05 (no irac) 
MAR05 
CB08 (no irac) 
CB08 



2.29Z 
2.29j 
2.29j 
2.29j 



2.29"* 



0.20: 

0.70l 
0.70l 
0.90^ 



1.00"* 



0.01] 

0.12^ 
0.10^ 
0.12"* 



0.15"* 



1-1-0.12 
-0.04 
1+0.19 

-0.04 
1+0.03 
-0.03 

^■^*-0.04 
2 94+0.03 
^•^*-0.03 

9 49+0.02 
^■^^-0.19 



+O.Vb 
1.85 
+0.75 
-1.05 
H.25 
-0.80 
1 7tr + 1.25 
^•'^-0.25 

^•'"^-1.00 

FO.OO 



+0.1U 
0.29 
+0.05 

0.24 
+0.20 
0.19 



1.9 

1.0 
0.4 
0.1 

1 2"'"^° 
4"*"° =' 



2.13: 



2.72"* 



o.o; 

11.9 
2.6 

1.1 

12.111°^ 
3.1 



0.19: 

0.58j 
O.60I 
0.78^ 



0.85"* 



0.1 
+0.5 

0.2 
+0.5 
0.1 



+ 1.25 
0.96 
+0.63 

-0.90 
+ 0.96 
0.33 



7-1-9.2 
-0.7 

,+1.5 
''-1.2 
,+1.9 
-0.6 

' 3 



+U.6b 
1.57 
+0.69 
-0.81 
+ 1.05 
0.64 



1.12 

1.06 
0.98 
0.93 
1.11 
1.15 



1256-142 



BC03 (no irac) 
BC03 

MAR05 (no irac) 
MAR05 
CB08 (no irac) 

CB08 



2.243 
2.24j 
2.24"* 



2.25 
1.75 
1.50 



0.30 
0.25 
0.20 



0.1 

0.2 
0.1 



2.76 

2.47" 
1.76 



0.7: 

1.2j 

0.6^ 



1.95 
1.50 
1.30 



'-'•^''-0.05 



1.00 



-0.20 - - -U.UM 

TTT7D n l9-l-i9-« 

-0.23 '^•-'-^-0.11 



0.20 
0.10 



T7U 

-0.23 
1.00 
0.75 
0.65 
0.56 
+0.00 
3.30 
0.50 
-0.23 
0.00 
0.2, 



+0.20 
0.19 
+0.02 
-0.09 



0"*"'* ° 
'^■'^-0.0 



0.4 
0.0 



+0.9 
0.4 
+0.1 
-0.0 



1 90"l"i-'*° 
^•^'^-0.35 



2.32 
1.61 



+0.2 
-0.7 
+0.3 
0.2 
+ 1.5 
0.6 
+1.5 
0.1 
+0.4 
0.3 
+0.0 
0.1 



+ 1.35 
0.71 
+0.31 
0.04 



1+0.78 
-0.60 
1+1.46 
-1.40 



0.1+°l 0.90 



444 7-^a4B.i 



1 50"^l-l'* 

1.0U_g 20 
1.0U_g 73 

+0.08 
-0.05 



,+0.82 
0.14 
,+0.89 
0.64 
+0.44 
0.43 
+0.04 
2.67 
+0.25 



1.19 
1.20 
1.27 
1.39 
1.10 
1.62 



1256-519 



BC03 (no irac) 
BC03 

MAR05 (no irac) 
MAR05 
CB08 (no irac) 
CB08 



1.857 
1.857 
1.857 
1.857 

1.857 
1.857 

,+0.19 
0.02 
+0.08 
-0.15 
+0.09 
0.14 
+0.09 
-0.15 
+0.08 
0.15 



0.30 

2.50 
0.60 
3.50 
0.30 
3.50 



o.yb 

-0.40 

0.45 

0.40 
0.60 
-0.30 
0.10 
0.70 
0.20 
-0.40 
+0.00 
0.30 



-0.11 
'^•^'^-0.39 
'^•J^>'-0.14 

'-'•^^'-'-0.09 
U.OU _Q 49 



3.2 
0.9 
2.2 
0.9 
3.0 
0.5 



4.64: 
3.51"* 



5.09 
4.24 
4.44 



+0.15 
-0.09 
+0.08 
0.11 
+0.10 
-0.09 
+0.02 
0.09 
+0.05 
0.09 
+0.03 
0.11 



TUTT 

0.5 
+0.4 
0.5 
+0.3 
-0.0 
+ 1.8 
0.0 
+0.5 
0.5 
+0.4 
0.1 



+0.34 
1.69 
+1.30 
1.20 
+0.43 
0.28 



+0.50 
0.22 
+0.18 
0.30 
+0.24 
0.13 
+0.85 
0.13 
+0.23 
0.28 

^•"^-0.02 



,+2.7 
3.0 
+252 
55.9 
+2.1 
4.8 

970 9+687.8 
^'O-^_276.0 
1 1+0.0 
^■-■■-l.l 



3.0 
55.9 

4.8; 



0.20: 

2.10 
0.46; 
2.86; 
0.21 



1+0.3 
-0.1 
1+0.3 
-0.1 
- + 1.4 
-0.5 



■--0.15 
3 00+° "*^ 



3+U.BU 
-0.32 
j+0.37 
-0.30 

0.801°;?° 

-0.09 



1.25 
2.55 
1.25 
2.33 
1.25 
3.07 



1256-1207 



BC03 (no irac) 
BC03 

MAR05 (no irac) 
MAR05 
CB08 (no irac) 
CB08 



1.66 
1.80 
1.80 
1.80 
1.80 



0.80 
0.80 
0.90 
0.90 
0.80 
1.00 



0.10 
0.12 
0.10 
0.10 
0.10 
0.12 



0.7; 

1.0 
0.0 
0.0 
0.6 
0.2 



1.40 
1.35 
0.78 
0.78 
1.24 



1.9"^ 

o.l^ 
o.l^ 

0.5"* 



0.70: 
0.681 



0.2 



2 ("19+0.25 
^■"^-0.18 
,+0.25 
-0.15 
7+0.00 
'-0.43 
3+0.04 
-0.02 
,+0.25 
-0.18 
,+0.25 
-0.16 



-0.00 



+0.01 
-0.00 
+ 0.06 
0.01 
+0.00 
-0.07 
+0.03 
0.04 
+0.01 
0.00 
+0.02 
0.09 



+0.2 
-0.0 
+0.4 
0.7 
+ 2.5 
-0.0 
+0.6 
0.0 
+0.2 
-0.0 
+ 1.7 
0.0 



+0.30 
-0.25 
+0.28 
0.52 
+1.15 
0.18 
+0.47 
0.01 
+0.18 
0.35 
+0.62 
0.30 



+ 1.3 
0.2 



0.80 
0.70 
0.88 



103.2 
-0.1 
+6.2 
-0.0 
+199 
-1.6 
+ 1.1 
-0.1 
1 4+98.3 
-'■•*-0.2 
1 9+0-6 



0.61 
0.19 
0.25 
+0.10 
0.23 



mwr 

-0.01 
+0.24 
0.09 
0.00 
-0.48 
0.07 
0.16 
+0.00 
0.01 
0.00 
0.51 



1.06 
0.99 
1.11 
1.18 
0.97 
0.95 



1256-1967 



BC03 (no irac) 
BC03 

MAR05 (no irac) 
MAR05 
CB08 (no irac) 
CB08 



2.02: 
2.27j 
1.99j 
2.02j 
2.02"* 



0.10: 
0.30^ 
0.60^ 
0.50 
0.10 
0.70 



0.10 
0.00 
0.55 
+0.10 
0.20 
+0.00 
0.02 
0.10 

-0.60 



0.01 
0.02 
0.08 
0.05 
0.01 
0.10 



1 

0.9 
0.0 
0.0 
1 

0.0 



2.43 
1.96 
1.85 
1.22 
2.43 
1.57" 



1.4 
0.0 

1-7: 

0.1 



0.09 
0.28 
0.52 
0.45 
0.09 
0.60 



-0.70 
+0.75 
-0.70 
h0.40 
-0.24 
h0.20 
-0.25 

1 oo+°-'''5 

l.UU_g 7Q 

in"'"° °° 

'-'■J^'-'_0.04 



5.0011^7° 

3+19.4 
-0.57 
3+0.00 
-19.9 
3+15.0 
-4.75 
1+0.01 
-0.03 



1+0.3 
"'-0.2 
+0.2 
-0.2 
+0.3 

0.3 
+0.2 

0.3 
+0.3 

0.2 
+0.0 

0.2 



+0.6V 
0.69 
+0.53 
-0.69 
+0.39 
0.58 
+0.15 
0.67 
+0.58 
0.69 
+0.04 
0.20 



+223.3 



+0.38 
-0.34 
+0.40 
-0.34 
0.19 
0.12 
0.10 
-0.11 
0.37 
0.34 
0.00 
0.01 



3.24 
4.15 
3.40 
3.23 
3.24 
4.24 



hdfsl-259 



BC03 (no irac) 

BC03 

MAR05 (no irac) 
MAR05 
CB08 (no irac) 
CB08 



2.249 
2.249 
2.249 
2.249 
2.249 
2.249 



1.00 
1.00 
0.30: 
0.30 



O.6O: 
20.0j 
5.00j 
0.04"* 



2.10 
2.10 
1.36 
1.38 
2.10 
1.10 



241.3 



-94.6 
,+140.0 
-93.0 

417 n+-*92.2 

*1'-U_246.3 

3+357.0 
'-280.1 
241.0t?j3^.5^ 
., + 158.5 
-261.1 



241.1^ 



542.8: 



289.6^ 



0.51 
0.51 
0.16 
0.15 
0.51 
0.06 



1.79 
1.62 
1.78 
1.63 
1.80 
1.83 
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TABLE 1 — Continued 



Name Model z age t A„ M* SFR < tsFR > < Xr > 

Gyr Gyr mag 10"Mq Mq yr'^ Gyr 

(1) (2^ (3) (4) (5) (6) (7) (8) (9) 

hdfsl-1849 BC03 (no irac) 2.30^^-^^ 0.401: JJ'^;; O^O^^r 1.5l^-? sITIJ^Y 0.32^^,11 Ta7 
BC03 2.30t°-°8 0.40t«'« COSl"-;? 1.61°^ 3.28lJ-i^ 35.5l|tl 0.32t°oit 1-34 
MAR05 (no irac) 2.30t° °! 0.40i»-5'^ Om+'^'H 0.9+li 1.96+111 1.7lf/ 0.35l° « 1.43 
MAR05 2.30t°°^ 0.40l°|° O.OSIOJ" O.qIJ'^ l-SetlH l-Vll^^'" 0.35^?,^? 1.33 
CB08(noirac) 2.301°;°^ OAOttll 0.081°^? l.sigj 2.961^;^,^ 32.0^11%' 0.32^^'^ 1.45 
CB08 2.26tlZ O-SOtgso 0-15lg;g? 0.9lg| 2.67l°;l^ H-SI^^ 0.65lg:| 1.39 

hdfs2-509 BC03 (no irac) 2.918 0.60+^'^^ 0.10+^-^^ 0.1+^'^ 4.28+6'jo' 14.41^-^ O.SO+U'Ij^, Oi 

BC03 2.918 0.50l°;J° 0.10l°;°2 0.4t°| 5.20tJ;|^ 45.7t|f;? 0.40tu;o5' 4.49 

MAR05 (no irac) 2.918 0.40lg;Jg 0.08t{];gg 0.4^°;° 4.53tg;gg 48.01^;^ 0.32t[J;of, 5.13 

MAR05 2.918 0.40lg'j° 0.08t°-°g 0.4t°-^ 4.5lto|/ 47.9t°j°g 0.32lj;''^« 4.69 

CB08(noirac) 2.918 0.60lg-°° O.lOig'gg O.llg'" 4.29l°-°° 14.3t°-° 0.50lg-°° 3.87 

^ CB08 2.918 O-eotgoo O-^O-o.oo ^-^-ofo ^■^'^-o-oo ^^-^^ao °-5°-.q.oo ^'^^ 

hdfs2-1099 BC03 (no irac) 2.72j:^-^^ CSOl^ i^ 0.12j:^'^^ 0.21^-^ 2.64^^-^1 3.7^^^ 0.68l^-J^ Tie 

BC03 2.5llOf,t 0.70t°-^o 0.12l0f8 0.7t°oi 3.09l°-g 10.li|-^ 0.58t1;ll 1.66 

MAR05 (no irac) 2.72l'^-'* 0.70tE;',V(! O.lOlEJf^ O.Ol?/; 2.04i0-g6 2.5±»-« 0.60tf,'« 1.44 

MAR05 2.55«;j° 0.701;',:- O.lOl-^.g O.ot^^ 1.68l°4^o 2.01;',:^ 0.60i;;:'« 1.54 

(..iJ()8 (110 ir;ic) 2.72_(, .;j 0.8()_ (, ().12_(, "-^-o 2 ^-O'lio .-,4 •lo_..2 1 ().o8_ (, (,<, 1.45 

hdfs2-2046 BC03 (no irac) 2.24lo JJ4 O.SOtlJ'fo 0.0210 01 O.SIq | I.SSIq 20 0.0+„o 0.28tiJ'oq ToO 

BC03 2.24l'^-E'| 0.40t;;-f(; 0.05l°f,i 0.5l'(J-| 1.27i;jf,i l.li?i 0.35i;;'j^ 0.97 

MAR05 (no irac) 2.24^0-;;^ O.SOlf/ftj 0.02l'^-;;5 O-^-ot l-26lo '1 0-0±o 0-27loo9 1-14 

MAR05 2.25ini 0.20l°-i° 0.011^^ 0.9l°i 1.43l0i? O.olg'g 0.19iro? 1.17 

CB08(noirac) 2.241001 0.40t0|0 O.OSIOO^ O.Stli 1.2(it°il ^■'^-H 0.35l:°V8 1-05 

^ CB08 2.24t°;°^ °-^°-o:io °-°^-o'o3 ^'"'^-"■li ^'"tO:^ ^■^^-o'?! ^'^^ 

ecdfs-4454 BC03 (no irac) 2.351 O.gol^'jo 0.201^;^^^ 0^2^^^!^ ^-^^^-o.w lO^T^o? 0.7ll^;^iJ Oi 

BC03 2.351 l-OOtl'll 0-30lo.io ^-^-ol 2.04l0;2? 32.9l2S;o 0.73lo;24 1.47 

MAR05 (no irac) 2.351 0.60tg jg 0.121°;°^ 0.2t° l l.OStlH 7.9ll_^i* 0.48tg;J| 1.45 

MAR05 2.351 0.60t°';i° 0-12lo'o2 0-2lo'2 I-O^lo'ol 7.9tJ-g 0.48t°-°| 1.34 

CB08 (no irac) 2.351 Cgol^'^o 0.20l:0-J0 0.2l:0-^ ^-'^f^ttfo 10.7lf,V O.nf^iil 1.57 

^ CB08 2.351 O-gO-ojo O-^O-ojs "-^-ol ^'^^-o'lfl ^O-^^jl °-^^-o:i5 ^'^^ 

ecdfs-4511 BC03 (no irac) 2.122 slxF^^^^ lO.Oli,^^ 1.2i;;-'j' Isl^^J^^^ 164.8l|?-| 1.57t|;'ii5 L98 

BC03 2.122 3.0015-00 20.0l08°° 1.210-0 4.0810^^7 168.3l«:2 l-53j:0'25 1.73 

MAR05 (no irac) 2.122 1-75^10, 20.0t0g00 1.2lo:2 2.76^0:97 i92.3t\l%^ 0.88lo:64 ^ gg 

MAR05 2.122 l-SOlj'^^ 5.00t^^g0 1.2+0-3 2.6010^^? 188.51^^*^1® O-^Sto'lg 1-85 

CB08 (no irac) 2.122 3.00lo:oo 20.0l0^00 ^ 2^0:0 4.07to;|0 167.3t«-9 1.53lo:i^ 2.01 

^ CB08 2.122 l-OOtg"" O-^^tojl l-llpl 2.37toi3 128.0^631 0-62lg;'?f 2.14 

ecdfs-4937 BC03 (no irac) 2.309 0.50t|;-oo 020^!^^™ LOl'/,-" S.OSl*;-"" ITiT^^^^Jl 0.34l|;-OjJ 1^28 

BC03 2.309 0.40tO;io 0.10l'^;0^ 0.7l'(j;g 2.471'^;^^ 58.3lJ03:6 o.30lO;'^| 1.30 

MAR05 (no irac) 2.309 0.20j;;5-i'; 0.05i5]-g^ 0.9tg-j 2.04t0-32 92.31^^^8® 0.15lg;g^ 1.25 

MAR05 2.309 0.lot°;ll 0-03t°i°ol ^-^-ol ^-^^-oli 252.9l?^9'*o 0.07lg:Jf 1.26 

CB08 (no irac) 2.309 0.50l0;0o 0.15t'^;0'^ 0.7t^;^ 2.54^0-^^ 79.8l|^;^ 0.36lg;?J 1.29 

^ CB08 2.309 O-SO-o'oo °-^°-o'po ^■'^-O-O ^•^^-6-io ^^"^-i^s ^'^^ 

ecdfs-5856 BC03 (no irac) S^SeT^J^^Il 0.50^0;^^ 0.08+^;^^ 0.7+^;| 1.93l^;i^^ ^^^^ 0.42lo-33 092 

BC03 2.55lo-J,|^ 0.60tfj!; 0.101J;:E;| 0.8+°ol 2.23tg;J| 7.51^^5^ 0.50tE;;^^ 1.05 

MAR05 (no irac) 2.56l'^-J| 0.60t0-^'(] O.OSl^j-^J^ 0.2tJ-^ 1.36tg-|^ l-2l?-2 0.52lO-J| 1.09 

MAR05 2.5610-00 0.30i0-« 0.03^0-07 o.8lO-| 1.69l0-g 0.31J0-9 o.27l°-33 1.O8 

CB08(noirac) 2.56l0-00 0.50lO-« 0.0810-0^ o.7lO-^ 1.9310-^1 e.lli^i^ 0.4210-26 0.91 

^ CB08 2.56lj;-0| 0.80l[',:| 0-12lg.g! O.llg.^ ^•^^-O-.i? '^■■^^^j "'^^-He ^-05 

ecdfs-6842 BC03 (no irac) 2.40lo-o{ O.SOl"-^;; o.25lo-3^ 1^2^^^! 3.741^-"^ 82^^^^ O^SS^^If^^ 2si 

BC03 2.40lo:o7 0.6010-^0 0.12lo:o8 0.810^1 2.32lo:|5 17,4+4.7 o.ASt^^ 2.81 

MAR05 (no irac) 2.40lo:o« 0.40l0-™ Cioio-jo Ll^O-S 2.44lo:68 QG.btlfo'' 0.30lO-J^ 2.60 

MAR05 2.40l0f7 o.20lo-jo 0.0410-01 l.slO-O 2.1510-03 44.4l0:/i 0.16lo-O0 2.60 

CB08 (no irac) 2.4010-?^; O.SOlo-^o 0.2510-^5 i,2l0-5 3.7310-06 81.61^03-^ o.58lo-i| 2.62 

CB08 2.43lg;g^ 0.70lg;gg 0.12lg;gg 0.2lgj 1.70lg;g 5.5lj:^ 0.58lji:;^ 3.52 

ecdfs-6956 BC03 (no irac) 2.037 1.501^;^^ 2(10^^^ l.Ol^;^ 1.43l^;:|i^ 118.8l^g;| 0.751^;!^^ 2^01 
BC03 2.037 l-25lj;^^ lO.Ol^o.o i.oio^a 1.25lg;|| 120.3l|g;^ 0.63i;;:|^ 1.97 
MAR05 (no irac) 2.037 0.60l0;^o 2.50li''7;; ^■'^-o.l 0-93lO;iI 168.8l^J;| 0.3ll0 ^f, 1.99 
MAR05 2.037 0-50lg;^',j 1.50l^^20 ^■'^-o.l 0-83l'o:o8 171.3l^g;f 0.26lg;',jg 1.89 
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TABLE 1 



Continued 



Name 
(1) 



Model 
(2) 



(3) 



age 
Gyr 
(4) 



Gyr 
(5) 



mag 
(6) 



M* 
10"Mq 
(7) 



SFR 
Mq yr- 
(8) 



< tSFR > 

Gyr 
(9) 



<Xr> 



CB08 (no irac) 
CB08 



2.037 
2.037 



1.50 
0.09 



-1.10 
0.01 
-0.03 



20.0 
0.03 



-0.55 
40+° ®5 

0.06l„„i 

10 ' 



+0.00 
19.2 
+0.02 
0.02 



1.0 
1.3 



^+U.UU 

-19 9 

2o.oi?r9 

^+18. 5 

-1.48 
,+ 19.8 
-0.10 
1+19.9 

"--(). f)f) 



+0.3 
0.1 
+0.1 
0.0 



1.43 
0.57; 



+U.4 
0.6 
+0.2 

0.3 
+0.6 
0.6 
+0.0 
0.0 
+0.4 
-0.7 
■ (1-1 
-0.8 



+0.35 
0.58 
+0.01 
0.09 



118.1 

118.7 



+U.84 
0.54 
+0.42 

0.38 
+0.29 
0.44 
+0.06 
0.11 
+0.68 
-0.(17 
4:! 
-O.l.'i 



+ 112.2 
-28.2 
+94.7 
103.9 



0.75 

0.06 



^316. 2 
-"-189.5 

154 q+509.9 
104-y_117.3 

+217.9 
147.6 
+533.3 

-199.4 



-0.55 
-0.00 
-0.01 



0.26 
+0.43 
0.16 

+0.50 
0.24 
+0.01 
0.01 

r,+0.74 

-().()() 



2.01 
2.22 



ecdfs-9822 



BC03 (no irac) 
BC03 

MAR05 (no irac) 
MAR05 

CB08 (no irac) 
('1308 



1.612 
1.612 
1.612 
1.612 

1.612 
1.012 



0.60^ 



20.0"' 



1.50: 
0.12^ 

0.01^ 



1.9 

2.0 

1-7; 

2.2 

1.9 
2.1 



0.96 
0.81 
0.75 
0.44 

0.49 
0.58 



196.0Z 



443.1 
15.0 
211.0 



0.30 

o.2o: 

0.26; 

0.04; 

0.05 
0.0() 



+0.09 
-0.05 
+0.08 
0.04 
+0.05 
-0.02 
+0.11 
0.03 
+0.13 
0.04 
+0.11 
-0.04 



0.1(1 
-0.10 
0.20 
0.00 
0.00 
-0.10 
0.00 
0.35 
0.10 
-0.10 
0.20 
0.00 



+0.03 
-0.04 
+0.08 
0.01 
+0.01 
-0.03 
+0.00 
0.07 
0.03 
-0.04 
+0.08 
0.01 



-FTTT- 

-0.5 
+0.0 
0.4 
+0.4 
0.1 
+ 1.1 
0.1 
+0.2 
-0.5 
+0.1 
0.1 



+0.12 
-0.20 
+0.16 
0.11 
+0.10 
0.06 
+0.20 
0.18 
0.08 
-0.21 
+0.17 
0.04 



4.5"^ 

+ 20.2 
13.0 
f 3.0 
-0.9 
+80.3 
2.4 
+3.8 
4.4 
+5.9 
0.0 



()..' 

-0.02 



0.0',- 
-0.00 
+ 0.14 

0.00 

0.00 
-0.08 
+0.00 

0.28 
+0.07 

0.08 
+0.12 

0.00 



1.30 
1.21 
1.30 
1.25 

1.32 
1.22 



ecdfs-11490 



BC03 (no irac) 
BC03 

MAR05 (no irac) 
MAR05 
CB08 (no irac) 

CB08 



2.37 
2.33 
2.33 
2.31 
2.33 
2.33 



0.30 
0.20 
0.30 
0.40 
0.30 
0.30 



0.05 
0.04 
0.04 
0.08 
0.05 
0.02 



0.5 
0.9 
0.1 
0.2 
0.5 
0.1 



0.70 
0.80 
0.54 
0.65 
0.67 
0.50 



16.9 
0.9 
6.9 

4.4; 

0.0 



0.25 
0.16 
0.26 
0.32 
0.25 
0.28 



0.30 
0.00 
-0.01 
+0.10 
0.22 
+0.00 
0.00 
+0.40 
0.30 
0.80 
-0.30 



^- '''-2.60 

o.oiinS 

1.50tl«4? 

,+0.00 

-0.00 



-I 9+0.0 
-'-■^-0.2 

1 o+O.l 
^••^-0.2 
-1 7+0.0 

i-'-o.o 

-,+0.0 

-0.2 
,+0.0 
^-0.2 



0.68ing 

„ oQ+0.07 

U.Ci9_Q QQ 

1+0.33 

-0.31 
-+0.60 
'-0.26 



266.3+[;-5 

-,+160.2 
-0.0 

388.7+"'^'^ 



0.14 
0.00 
-0.00 
+0.06 
0.10 
+0.00 
0.00 
+0.19 
0.13 
+0.41 
-0.14 



2.12 
2.43 
2.06 
2.38 
2.12 
2.31 



ecdfs-12514 



BC03 (no irac) 
BC03 

MAR05 (no irac) 
MAR05 
CB08 (no irac) 
CB08 



2.024 
2.024 
2.024 
2.024 
2.024 
2.024 



0.60 
0.04 

o.3o: 

0.02 

0.60 
0.70 



145.2^ 



0.01"' 



3.50 
20.0 



0.70 
+0.70 
0.50 
+0.70 
0.17 
+0.20 
0.15 
+ 1.20 
0.70 
0.91 
-0.02 



+ 16.5 

-3.30 
+0.00 
-19.2 



1.2"' 

l.2^ 



i.4i: 

1.55j 



+ 19.6 
0.37 
+0.40 
0.22 
+19.9 
0.09 



"'-0.8 
,+0.4 
^-0.5 
O+0.7 
^-1.0 
O+0.2 
'-0.9 
+0.8 
0.7 
+0.0 
-1.3 



+ 1.02 



191.8 
+ 124.8 
0.0 

+75.4 
115.7 

271.lt||j|_ 



660.3 

268.1 



0.31 
0.03 
0.15 
0.01 

0.30 
0.35 



+46,-). 1 
74.8 
+86.9 
38.1 
zio. u_]^9g 3 
91 9 9 + 131.8 

1+80.4 
0.8 



HU.bl 
-0.38 
-+0.57 
-0.34 
,+0.27 
-0.10 

hO.49 
-0.45 
+ 0.72 
-0.01 



2.10 
3.36 
2.12 
3.15 

2.09 
2.79 



ecdfs- 13532 



BC03 (no irac) 
BC03 

MAR05 (no irac) 

MAR05 
CB08 (no irac) 
CB08 



2.336 
2.336 
2.336 
2.336 
2.336 
2.336 



0.80 
0.80: 
0.20 

0.20; 

0.80 
0.09 



0.40 
0.30 
O.IO 



10+0-20 
^•1^-0.09 



0.40 
0.01 



+0.U2 



1-0.02 
-0.00 
3+0.00 
-0.00 
,+0.01 
-0.00 
,+0.00 
==-0.00 

03+"'"^ 

^■^''-0.00 
FO.OO 
-0.00 



+ 19.6 
0.37 
+0.24 
-0.00 



1, 
1, 
1, 
1, 

1.4 
1.7 



1.75 



1.64"^ 



-0.71 
,+0.27 
-0.39 
1 i4+°-52 
^■^^-0.53 
1 14+0.20 



1.75 
0.96 



01+"-"^ 
u-'JJ^-o.oo 

1+0.00 
-0.00 
1+0.29 
-0.00 
1+0.01 
'--0.00 
I +0.02 
-0.00 

''•'^^-0.00 



,+0.0 
-0.0 
,+0.0 
-0.0 
,+0.0 
-0.0 
,+0.0 
^-0.0 
1+0.0 
-0.0 
1 1+0.0 
-0.0 



-0.46 
+0.82 
0.71 
+0.33 
-0.16 



87.9 
52.6 



1.4-' 



0.52: 
0.55^! 
0.131 



0.52 
0.08 



o.68in? 

0.68+-°° 



"396!oI|^^ 
396.21?-^ 

865.4t«f/ 
396.3t;!;;; 



0.00 
0.00 
+0.00 
0.00 
+0.00 
0.00 
+0.00 
0.00 
0.00 
0.00 
0.00 
-0.00 



0.97 
0.88 
0.96 

0.87 
0.98 
1.02 



ecdfs-16671 



BC03 (no irac) 
BC03 

MAR05 (no irac) 
MAR05 
CB08 (no irac) 
CB08 



2.61 



-0.01 
,+0.00 
^-0.01 
,+0.04 
^-0.02 
,+0.04 
-0.02 
, +0.04 
-0.01 
2 61+° °° 



2.6i: 

2.63j 
2.63j 
2.61^ 



0.03: 

0.031 

0.021 

0.02^ 



0.03^ 



0.01: 
O.Oll 
0.01^ 

0.01^ 



1.1- ' 

l.l^ 

1.2- * 
1.2-* 
l.l^ 



irw 

-0.40 
0.40 
0.30 
0.20 
0.30 
+0.10 
0.14 
0.30 
-0.30 
0.00 
0.01 



+0.00 
19.6 
+ 19.3 
0.50 
+0.00 
19.9 
+3.35 
0.62 
+ 15.0 
-4.60 
+0.01 
0.00 



q+O.i 

■•-'-0.2 
9+0.0 
•^-0.1 
4+0-2 
■^-0.2 
c+O.O 
■>'-0.2 

q+0.1 

■•-'-0.2 

q+0.1 

■•^-0.0 



-0.00 
,+0.05 
-0.00 
-0.0 

■--0.01 
j+0.06 
-0.03 

n fis+0-00 



0.51: 

0.5lj 

0.68^ 



0.02 
0.02 

0.01 
0.01 
0.02 
0.02 



+0.38 
0.55 
+0.32 
0.42 
+0.37 
0.67 
+0.30 
0.44 
+0.51 
-0.38 
+0.00 
0.04 



356.8l'^|4=:i 
259.8111:1 

51 5 q+266.7 



-0.12 
-0.15 
0.15 
0.15 
+0.10 
0.14 
+0.06 
0.06 
+0.16 
0.12 
+0.00 
0.00 



4.71 
4.46 
4.90 
4.64 
4.70 
4.45 



cdfs-6202 



BC03 (no irac) 
BC03 

MAR05 (no irac) 
MAR05 
CB08 (no irac) 
CB08 



2.225 
2.225 
2.225 
2.225 
2.225 
2.225 



0.80 
0.60 
0.40 
0.20 
0.70 
0.10 



20.0 
0.65 
20.0 
0.65 
5.00 
0.03 



2.32 
2.04 
1.74 
1.30 
2.15 
1.14 



649.3 
357.8 
I68.O: 



+27.1 

280.2 
+93.7 
-165.3 
-121.4 
-0.0 



0.40 
0.34 
0.20 
0.10 
0.35 
0.07; 



1.30 
1.87 
1.35 
1.97 
1.31 
2.39 



TABLE 2 

CoMrWRT^ON OF PHOTOMETRIC ReD^IIIFT^ 



Data Set 


Ofltsct 


Scatter 




(^phot-2spec)/(l H 


h z) (7/(1 + 2) 


All U^K 


0.17 


0.12 


All U^4.5/im 


0.18 


0.12 


All U^8.0/im 


0.14 


0.11 


Deep U^K 


0.00 


0.05 


Deep U^4.5^m 


0.01 


0.05 


Deep U— +8.0^m 


0.02 


0.05 


Wide U^K 


0.46 


0.19 


Wide U^4.5/im 


0.46 


0.18 


Wide U-»8.0/nm 


0.35 


0.17 
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TABLE 3 

Systematic Offsets in SED Parameters Compared to 
U— ►S^iM+NIRsPEC Fits 



Data 


Log(M0) 


T 

Log(Gyr) 


{*)SFR 
Log(Gyr) 


mag 


SFR 
Log(M0 yr-i) 


U-»K 

U->z'+NIRspec 


0.111 ± 0.045 
0.022 ± 0.020 
-0.012 ± 0.019 


0.091 ± 0.138 
0.039 ± 0.114 
-0.055 ± 0.154 


-0.009 ± 0.078 
0.002 ± 0.054 
0.015 ± 0.062 


0.070 ± 0.137 
0.032 ± 0.094 
0.026 ± 0.055 


-0.003 ± 0.187 
-0.024 ± 0.118 
-0.080 ± 0.076 



Note. — Quoted values are the mean of Log(Data/U— >8/im-|-NIRspec), except for A„ which is mean (Data - U— >8/^m-|-NIRspec). 
Errors are standard errors of the mean. 



TABLE 4 

Effect of Stellar Synthesis Code on SED Parameters 



Parameter 




MAR05/BC03 


All 




CB08/BC03 




EL 


NEL 


EL 


NEL 


All 


T 


0.66 (0.84) 


0.83 (0.84) 


0.75 (0.84) 


0.66 (111.) 


1.25 (1.48) 


1.00 (59.8) 


(*>SFR 


0.42 (0.64) 


1.03 (1.00) 


0.65 (0.81) 


0.96 (1.39) 


1.46 (1.53) 


1.24 (1.45) 




0.10 (-0.0) 


-0.4 (-0.3) 


0.00 (-0.1) 


-0.2 (-0.1) 


-0.6 (-0.5) 


-0.4 (-0.3) 




0.63 (0.65) 


0.69 (0.69) 


0.63 (0.67) 


0.70 (0.76) 


0.77 (0.77) 


0.74 (0.76) 


SFR 


1.58 (1.79) 


0.33 (0.69) 


1.12 (1.27) 


0.64 (1.04) 


0.31 (1.64) 


0.32 (1.33) 



Note. — Quoted values are the median ratio, with the mean ratio in brackets 



TABLE 5 

Effect of Dust Law on SED Parameters 



T 

(*>SFR 
Au 

SFR 



Parameter BC03+SMC/BC03+Calz 



EL 



NEL 



All 



BC03-l-LMC/BC03+Calz 
EL 



NEL 



All 



BC03-l-MW/BC03+Calz 
EL 



NEL 



All 



2.00 (16.7) 
2.81 (6.62) 
-0.6 (-0.6) 
1.20 (1.37) 
0.55 (3.82) 



3.33 (131.) 

2.35 (5.56) 
-0.1 (-0.3) 
1.07 (1.34) 

2.36 (-58.5) 



2.50 (70.7) 
2.80 (6.12) 
-0.4 (-0.5) 
1.15 (1.36) 
1.00 (29.5) 



1.00 (2.73) 
1.11 (1.07) 
-0.2 (-0.1) 
0.91 (0.80) 
0.77 (210.) 



1.25 (1.33) 
1.00 (1.09) 
-0.1 (-0.1) 
0.94 (0.91) 
1.28 (2.75) 



1.00 (2.07) 
1.09 (1.08) 
-0.2 (-0.1) 
0.92 (0.85) 
1.00 (112.) 



1.00 (1.60) 
1.00 (0.93) 
-0.1 (-0.1) 
0.83 (0.73) 
0.73 (208.) 



1.00 (1.14) 
1.00 (1.06) 
-0.1 (-0.1) 
0.93 (0.91) 
0.93 (2.35) 



1.00 (1.38) 
1.00 (0.99) 
-0.1 (-0.1) 
0.90 (0.82) 
0.83 (111.) 



Note. — Quoted values are the median ratio, with the mean ratio in brackets 



TABLE 6 

Effect of Metallicity on SED Parameters 



Parameter 


BC03.m42/BC03.m62 


BC03.m72 /BC03.m62 




EL NEL All 


EL NEL All 



T 0.92 (0.84) 0.39 (0.97) 0.50 (0.90) 0.53 (28.3) 1.25 (126.) 0.79 (74.3) 

(*>SFR 0.47 (0.68) 1.17 (1.36) 0.94 (1.00) 0.86 (1.75) 1.17 (1.61) 1.03 (1.68) 

A^ 0.39 (0.43) 0.19 (0.06) 0.20 (0.26) -0.3 (-0.3) -0.4 (-0.4) -0.3 (-0.3) 

Mstar 0.74 (0.71) 1.03 (1.01) 0.91 (0.85) 0.83 (0.95) 0.81 (0.91) 0.82 (0.93) 

SFR 1.04 (180.) 0.39 (1.43) 0.65 (96.4) 0.59 (1.60) 0.67 (1.53) 0.59 (1.57) 



Note. — 



Quoted values are the median ratio, with the mean ratio in brackets 



